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A PROBABILISTIC  ANALYSIS  OF  EMBANKMENT  STABILITY  PROBLEMS 


by 

LAWRENCE  WILLIAM  GILBERT 


Submitted  to  the  Department  of  Civil  Engineering 
on  (August,  1974)  in  partial  fulfillment  of  the 
requirements  for  the  degree  of  Master  of  Science 
in  Civil  Engineering. 


In  this  thesis,  a probabilistic  model  is  developed 
to  predict  the  reliability  of  an  embankment  constructed 
on  soft  saturated  clay.  The  model  is  based  on  a circular 
arc  method  of  analysis,  supplemented  with  a measure  of  the 
uncertainty  in  the  resisting  and  in  the  overturning  moments. 
The  uncertainty  in  the  overturning  moment  was  considered 
negligible  in  this  thesis.  The  uncertainty  in  the  resis- 
ting moment  was  considered  due  to  the  uncertainties  of 
bias,  random  testing  error,  and  inherent  soil  variability. 

Two  case  studies  were  analyzed  in  this  thesis  by 
both  the  conventional  method  of  analysis  and  the  proba- 
bility model.  The  results  indicate  that  the  uncertainties 
in  bias  correction  factors  are  the  dominant  sources  for 
both  field  vane  testing  and  unconfined  compression  testing. 

The  basic  probability  model  is  then  extended  to  in- 
clude the  effect  of  embankment  length  on  the  computed 
failure  probability.  Two  approaches  are  taken.  The  first 
is  a direct  extension  of  the  basic  model,  considering  the 
actual  embankment  length  as  a multiple  of  the  minimum  em- 
bankment length  required  to  satisfy  the  assumption  of 
"plane  strain."  The  second  approach  is  a three  dimensional 
probability  model  developed  from  a first  passage  failure 
criter i on . 
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Introduction 


Probabilistic  methods  only  recently  have  been  intro- 
duced into  the  field  of  geotechnical  engineering.  There 
has  been  some  skepticism  as  to  the  worth  of  probability 
models  in  the  field.  The  two  main  reasons  for  the 
skepticism  are: 

1.  Measurements  of  soil  properties  such  as  undrained 
shear  strength  (S^)  can  be  biased.  The  magnitude 
of  bias  has  to  be  evaluated  if  the  results  from 
probability  models  are  to  have  meaning. 

2.  The  scatter  of  soil  property  measurement  is  not 
purely  random.  The  geologic  history  of  a soil 
deposit  can  result  in  zones  of  soil  having  a 
strength  difference  from  the  average  for  the 
layer.  The  problem  of  test  bias  for  the  most 
common  strength  tests  (e.g.  the  field  vane  and 
the  unconfined  compression)  and  the  inherent 
variability  of  soil  properties  are  treated  in  a 
rational  manner. 

Models  developed  from  probability  theory  are  subject 
to  limitations,  as  are  all  engineering  models.  Proba- 
bility models  do  not  replace  "engineering  judgement". 

Rather,  they  provide  frameworks  that  allow  the  engineer 
to  exercise  this  judgement  in  a systematic  manner. 
Probabilistic  models  can  assist  one  to  evaluate  both  the 
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limitations  and  consequences  of  design. 

This  thesis  presents  a probability  model  for 
analyzing  the  undrained  stability  of  earth  embankments 

constructed  on  soft,  saturated  clay.  It  consists  of  five 
chapters . 

Chapter  One  presents  the  reasons  why  the  conventional 
safety  factor  does  not  indicate  the  true  reliability  of 
an  embankment.  It  also  presents  in  general  terms  the 
advantage  of  expressing  the  safety  of  an  embankment  by  its 
probability  of  failure  rather  than  its  conventional  safety 
factor . 


Chapter  Two  introduces  the  basic  concepts  and  tools  that 
will  be  applied  throughout  this  report.  Then,  a two- 
dimensional  probability  model  is  developed  for  an  embank- 
ment on  a homogeneous  foundation.  This  is  the  case  of  an 
embankment  having  an  assumed  failure  plane  within  a single 
layer  of  uniform  strength,  of  constant  inherent  varia- 
l)ility,  and  uniformly  tested.  Next,  the  model  is  extended 
to  include  layered  soils.  The  chapter  also  includes  a 
presentation  of  the  methodology  for  determining  the  input 
for  the  model. 

Cliapter  Three  analyzes  case'  studies  with  the  probability 
moilel.  It  as  the  intent  of  this  cliapter  to  clarify  the 
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model  by  applying  it  to  real  problems  on  which  typical 
soil  data  are  available.  Differences  between  the  conven- 
tional safety  factor  and  the  probability  of  failure  also 
are  illustrated  in  this  chapter. 

Chapter  Four  contains  an  extention  of  the  general  two- 
dimensional  model.  In  this  chapter  the  influence  of 
embankment  length  on  the  probability  of  failure  of  an 
embankment  are  evaluated.  The  examples  in  chapter  three 
are  used  to  illustrate  this  extension. 

Chapter  Five  presents  the  conclusions  from  the  study 
and  makes  recommendations  for  future  research. 


( 

I 

! 
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C liapter  1 Sji f e ty  F a cto r s and  F ailure  Probabilitie s 
1.1  Genera  I Faced  with  the  responsibility  of  designing 
an  earth  embankment  against  a shear  failure,  an  engineer 
practicing  the  current  design  approach  begins  by  selecting 
the  appropriate  method  of  stability  analysis.  Since  the 
probability  model  developed  in  this  thesis  is  for  embank- 
ments on  soft  saturated  clays,  stability  during  and 
immediately  after  construction  is  most  critical.  As  an 
embankment  is  constructed,  positive  pore  pressures  are 
induced  in  the  soft  foundation  clays.  As  the  excess  pore 
pressure  dissipates,  the  strength  of  the  foundation 
increases.  This  means  that  the  safety  of  the  embankment 
improves  with  time. 

The  appropriate  method  of  stability  analysis  for 
the  above  case  is  the  total  stress  analysis  (TSA) . The 
strength  parameter  required  for  a total  stress  analysis 
IS  the  undrained  shear  strength,  Su. 

Once  the  method  of  analysis  is  ascertained,  the 
engineer  seeks  information  on  the  soil  typos  and  layering 
geometry  under  the  proposed  embankment.  This  is 
accomplished  with  sub  surface  exploration  (and  soil 
class! f ication  of  recovered  samples)  at  the  site  of  the 
proposed  embankment.  Hopefully,  the  site  exploration  is 
<ii<lcd  by  a knowledge  of  the  geology  of  the  area. 

Next,  the  c'nginoer  tries  to  measure  the  undrained 
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strength  with  depth  and  embankment  length  along  the 
proposed  embankment.  The  undrained  strength  can  be 
measured  with  either  field  tests  or  laboratory  tests  on 
"undisturbed"  samples.  Current  engineering  practice 
usually  uses  either  field  vane  tests  or  unconfined 
compression  (or  triaxial  unconsolidated-undrained)  tests 
on  representative  "undisturbed"  samples. 

With  the  proposed  embankment  geometry,  a series  of 

total  stress  stability  analyses  are  performed  to  locate 

the  most  critical  failures  plane.  The  most  critical 

failure  plane  is  that  which  produces  the  minimum  safety 
( * ) 

factor.  The  design  engineer  would  compare  this  mini- 

mum safety  factor  to  a specified  "design"  safety  factor 
and  adjust  the  conf igurations  of  the  section,  if 
necessary,  until  the  safety  factor  of  the  section  is 
equal  to  or  slightly  greater  than  the  "design"  safety 
factor . 

Historically,  ranges  of  design  safety  factors  have 
been  used  to  account  for  unknowns  associated  with 
different  types  of  foundation  designs  (Meyerhof,  1970). 
These  design  safety  factors  reflect  values  that  have  led 

(*)  The  assumption  of  a circular  failure  mechanism  is 
made'  in  this  the'sis.  Based  on  this  assumption  the  safety 
factor  is  generally  defined  as  equal  to  the  resisting 
motm'nt  divided  by  the  overturning  moment,  i.e. 

FS  = MU 
Wo 
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to  satisfactory  performance  based  on  both  stability  and 
deformation.  The  magnitude  of  the  safety  factor  in  some 
way  reflects  uncertainties  due  to  simplified  soil  and 
embankment  configurations,  simplified  stability  analyses 
and  loads  and  resistances  generated  by  a given  structure. 

A set  of  design  safety . factors  was  proposed  by 
Torzaghi  and  Peck  (1967)  and  is  given  in  Tabic  1-1.  The 
lower  and  upper  limits  reflect  those  that  should  be  used 
for  maximum  and  average  design  loads,  respectively.  Also 
shown  in  Table  1-1  are  the  failure  rates  of  the  three 
classes  of  structuresas  estimated  by  Meyerhof  (1970). 


Class 

(T&P,  1967) 
Safety 
Factor 

(Meyerhof ,1970) 
Approximate 
Failure  Rates 
(per  1000) 

Earthworks 

1.3  to  1.5 

5 

Earth  Retaining  Structures 

1.5  to  2.0 

1 

Foundations 

2.0  to  3.0 

0. 5 

Design  Safety  Factors  and  Approximate  Failure  Rates 


TABLF  1 - 1 

1.2  Limitations  of  Current  Design  Approach 

The  current  design  approach  for  analyzing  the 
st.ibiiity  of  earth  emb.inkmcnts  on  soft  saturated  clays 
involvivs  three  basic  st«-ps.  They  are: 

1.  Solec't  an  appropriate  mi'thod  of  total  stress 
analysis,  the  seli'cted  method  usually  assumes 


-16- 


either  a circular  arc  or  a wedge  type  failure 
surface.  These  two  types  of  failure  mechanisms 
have  different  definitions  of  safety  factor. 

2.  Estimate  the  input  parameters  for  the  selected 
method  of  analysis.  The  undrained  shear  strength 
is  needed  to  calculate  the  resistance  to  failure. 
The  total  unit  weight  of  the  embankment  and  the 
foundation  geometry  are  needed  to  calculate  the 
overturning  moment. 

3.  Subdivide  the  embankment  length  into  sections  of 
similar  subsurface  conditions.  A separate 
design  is  made  for  each  section. 

There  are  four  basic  limitations  to  applying  the 
current  design  approach.  They  are; 

1.  Uncertainty  in  the  method  of  analysis.  There  is 
an  uncertainty  as  to  the  true  meaning  of  safety 
factors  obtained  from  conventional  methods  of 
analysis.  This  uncertainty  is  due  to  inconsis- 
tencies in  the  definition  of  safety  factor  and 
to  simplifying  assumptions  made  in  developing 
and  applying  conventional  methods  of  analyses. 

2.  Uncertainty  in  the  Resistance.  This  uncertainty 
is  due  to  two  sources.  First,  uncertainty  is 
introduced  in  tryin<)  to  mt'asure  the  true  in  situ 
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strength  of  the  foundation  clays.  Uncertainty 
due  to  bias  error,  random  testing  error,  and 
inherent  soil  variability  obscure  the  true  or  in 
situ  undrained  strength.  The  second  source  of 
uncertainty  lies  in  what  value  of  strength  is 
imput  into  the  analysis.  Some  engineers  use  the 
mean  value,  while  others  use  rc'duced  means. 

There  is  also  uncertainty  as  to  whether  the 
strength  of  fill  and  the  fill  crust  strength 
should  bo  considered. 

3.  Uncertainty  in  the  Load.  This  uncertainty  is 
due  to  uncertainties  in  measuring  the  total  unit 
weight  of  the  embankment  material.  Assuming 
that  the  embankment  is  constructed  on  level 
ground,  the  unit  weight  of  the  foundation  clays 
does  not  contribute  to  the  net  load.  Uncertainty 
in  tile  load  is  also  contributed  by  uncertainty  in 
the  embankment  conf iguration  and  in  externally 
applied  loads. 

4.  Uncertainty  due  to  embankment  length.  The  safety 
factor  of  an  embankment  is  generally  independent 
of  its  length.  If  two  embankments  of  different 
lengths  are  constructed  to  the  same  safety 
factor  on  the  same  foundation,  the  longc'r  embank- 
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ment  will  have  the  higher  failure  probability. 

These  limitations  will  be  explained  in  greater  detail 
in  the  following  section. 

Determining  the  safety  factor  of  an  earth  embankment 
requires  selecting  a method  of  analysis,  estimating  the 
loads  and  resistances,  and  subdividing  the  length  of 
embankment  into  sections  of  similar  subsurface  conditions. 
Uncertainties  are  associated  with  each  of  these  stops.  It 
IS  a quantitative  evaluation  of  these  uncertainties, 
together  with  the  embankment  safety  factor,  that  defines 
the  probability  of  failure  of  an  embankment.  The 
uncertainties  in  these  steps  will  be  treated  separately. 

1.3  Uncertainty  in  the  Method  of  Analysis. 

1.3.1  Safety  Factor  Definition 

Uncertainty  as  to  the  meaning  of  a safety  factor  is 
introduced  with  the  definitions  of  safety  factor.  Con- 
sider first  an  example  of  an  embankment  investigated  by 
the  wedge  method  of  analysis.  Using  the  notation  shown  in 


Figure 

1-la,  the  safety  factor 

against  si 

iding  is 

genera 

lly  defined  as: 

Ra  + Rb 

+ Rp 

ba- 

Dp 

(1.3-1) 

in  whu 

ch  D = driving 

force;  R 

- resisting 

force;  ft.  refers 

to  the 

act  1 ve  wf'dge  ; 

b refers 

to  the  cent 

ral  block;  P 

re  f ers 

to  the  passive 

wedge . 
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The  safety  factor  for  the  wedqe  method  of  analysis 


J 

4 


has  been  defined  as: 


PS  _ Ra  -*■  Rb  + Rp  + Dp 
Da  ’ 


(1.3-2) 


Whether  tiie  term,  Dp,  is  placed  in  the  numerator  or  in 


the  denominator  of  the  equation  of  the  safety  factor  will 
make?  a difference  in  the  calculated  safety  factor  of 
an  {'inbankment.  The  calculated  safety  factor  for  the  two 
equations  are  the  same  only  when  the  safety  factor  is 
equal  to  unity. 


A similar  inconsistency  exists  for  the  circular  arc 
method  of  analysis.  Considering  the  same  embankment  and 
the  notation  shown  on  Figure  1-lb,  the  safety  factor  is 
generally  defined  as; 


FS 


(1.3-3) 


where  S = shear  strength;  L = failure  arc  length;  r = 
radius  of  failure  arc;  a^^  and  - moment  arms;  and 
= weights  of  masses  1 and  2.  However,  the  safety 
factor  for  the  circular  arc  method  of  analysis  was 
initially  defined  by  Terzaghi  and  Pock  (1948)  as: 


r . e _ S L r + Wo  a o 

____  (1.3-4) 

Again,  the  location  of  the  term,  a^,  will  result 
in  different  values  of  safety  factor  for  safety  factors 
greater  than  unity.  Also,  there  is  no  reason  why  the 
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safety  factors  bash'd  on  a circular  arc  analysis  should 
*''jual  t.o  the  safety  factors  of  a slidinq  wedqe 
analysis  (Wricjht  and  Duncan,  1972)^  Thus,  an  embankment 
can  liavf'  at  least  four  different  safety  factors  dopendinq 
on  the  selected  method  of  analysis  and  on  the  definition 
of  safi'ty  factor  for  that  method  of  analysis. 

This  inconsistency  in  the  definition  of  safety 
l.ictor  can  be  minimized  in  a probability  model  based  on  a 
circular  <irc  method  of  analysis.  The  probability  model 
beqins  with  the  definition  of  safety  margin,  M,  or  the 
net  resisting  moment. 

M = safety  margin  = Mr  - Mo  (1.3-5) 

M = (SLr  + a^)  - a^  = 

M = SLr  - (W^  a^  - a^)  (1.3-6) 

Whether  the  term  W2  a^  is  considered  as  a negative  over- 
turnini)  moment  or  a positive  resisting  moment  docs  not 
ch.inge  the  safety  margin. 

1.1.2  .S  impl  1 fy  1 nj  Model  Assumptions 

Simplifying  assumptions  made  both  in  developing  and 
in  applying  stability  analyses  result  in  methods  of 
limited  reliability.  Some  of  these  assumptions  are: 

1.  Approximations  of  the  embankment  configuration, 
foundation  layering,  and  failures  surface 
shape  are  all  sources  of  uncertainty. 
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2.  The  behavior  of  the  embankment  may  not  agree 
with  the  assumptions  of  the  model  (plastic 
deformations  rather  than  rigid  body  motion, 
anisotropic  versus  isotropic  strength,  and 
strain  softening  material). 

3.  The  conventional  methods  of  analysis  assume 
plane  strain  conditions.  Many  observed  failures 
are  bowl-shaped,  far  from  plane  strain. 

1 . 4 Uncertainty  in  the  Resistance 
1.4.1  Measurement  of  In  Situ  Strength 

Since  this  report  is  concerned  with  the  undrained 
stability  of  an  embankment  constructed  on  soft  saturated 
clays,  the  soil  property  needed  for  calculating  resistance 
is  the  undrained  shear  strength.  Selecting  the  in  situ 
strength  is  usually  a major  source  of  uncertainty  in 
stability  analyses.  The  uncertainty  in  measuring  strength 


will  be  illustrated  with  the  field  vane  test  as  an 
example.  The  basic  problems  with  other  lab  and  field 


tests  are  similar  to  those  of  the  field  vane. 

The  first  problem  with  the  field  vane  test  is  bias. 

A lest  IS  biased  when  it  consistently  overpredicts  or 
undi ■ r pred i ct s the  value  of  in  situ  strength,  n>gardless 
of  the  number  of  tests  performed.  There  are  six 
d i sadvantaejes  of  the  field  vane  test  that  lead  to  biased 


-22- 


K 


moasuroment  (Ladd,  1971).  They  are: 

1.  The  stress  system  during  shear  is  unli)ce  any 
typical  stress  system  applied  to  soils.  Since 
the  undrained  strength  of  a soil  is  dependent 
upon  the  stress  system  at  failure,  the  field 
v.im-  can  only  be  a semi-emperical  measurement 
ot  strength. 

2.  .Sample  disturbance  during  testing  leads  to  an 
underestimation  of  the  in  situ  shear  strength. 

For  a series  of  tests  on  quiclt  clay,  field  vanes 
left  in  the  ground  for  one  day  after  penetration 
gave  substantially  higher  strengths  than  those 
sheared  immediately  after  insertion.  (Flaate, 
1966).  This  is  a measure  of  the  disturbance  due 
to  the  method  of  testing. 

3.  The  presence  of  sand  lenses,  pieces  of  wood, 

or  shells  can  increase  the  scatter  of  field  vane 
tests  and  can  increase  the  mean  of  the  tests 
significantly.  The  field  vane  test  is  best 
suited  for  " homogi'neous " soil  deposits. 

4.  Since  the  test  is  performed  in  the  field,  close 
supervision  by  the  engineer  is  essential.  This 
Will  assure  that  the  tests  are  performed  in  a 
standardized  manner  regardless  of  the  drill  crew. 
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0.  The  field  vane  test  is  not  applicable  to  very 
stiff  clays . 

6.  Partial  drainage  may  occur  in  performing  this 
field  vane  test  thus  altering  the  measured 
strength . 

These  sources  of  bias  can  be  considered  in  part  for 
by  tleveloping  a correction  factor.  By  correlating  field 
vane  measurc!ments  to  "in  situ"  strengths  back  calculated 
from  embankment  failure,  a correction  factor  can  be 
obtained  by  relating  the  "in  situ"  to  the  measured  field 
vane  strength  (Bjerrum,  1972).  The  product  of  the 
measured  strength  and  the  correction  factor  yields  the 
best  estimate  of  strength. 

This  allowance  for  test  bias  is  essential  to  bot.h 
the  conventional  method  of  analysis  and  the  probability 
model.  llowc.-v<?r,  the  conventional  method  does  not 
explicitly  allow  for  the  uncertainty  associated  with  the 
development  of  a correction  factor.  The  reliability  of 
such  strength  corrections  can  have  a significant  effect 
on  the  true  probability  of  failure  of  an  earth  embankment. 

The  second  problem  with  the  field  vane  test  is 
random  testing  error.  Even  if  bias  is  completely 
eliminated  the  individual  test  results  will  be  scattered 
about  the  mean  stren<jth  due  to  testing  errors.  Increasing 


-24- 


I 


scatter  implies  an  increase  in  the  probability  of  failure 
of  an  embankment.  If  a test  had  no  bias  and  no  random 
error  it  would  be  a perfect  test,  exactly  measuring  the 
in  situ  strength  at  a point.  For  this  case  there  would 
bo  no  uncertainty  in  measured  strength.  In  general,  the 
variability  of  test  results  from  field  vanes  is  high. 

This  uncertainty  can  have  a significant  effect  on 
the  probability  of  failure  of  an  embankment,  but  can  have 
no  direct  effect  on  safety  factor  determinations. 

However,  some  design  engineers  may  use  lower  values  of 
strength  than  normal  to  account  for  increased  scatter. 

The  third  problem  is  soil  variability  - how  in  situ 
strength  varies  with  dimensions.  This  problem  is  due  not 
to  the  type  of  test  performed,  but  to  the  soil  itself. 

If  the  field  vane  test  were  a perfect  test,  having  no 
bi.is  nor  random  error,  the  scatter  of  test  results  would 
I reflect  inherent  soil  variability.  As  inherent  soil 

I variability  increases  the  [irobability  of  failure 

I 

generally  increases,  but  the  conventional  safety  factor 
can  be  inaffected. 

i 

“ The  three  problems  with  field  vane  testing  that 

> introduce  uncertainty  into  stability  calculations  are 

c'ach  illustrated  on  Figure  1-2  and  1-3. 

Fiijure  1-2  demonstralc's  the  fjualitative  importance 
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of  random  testing  error  and  inherent  soil  variability. 

The  first  case  presents  the  determination  of  strength  for 
the  same  subsoil  by  two  different  unbiased  testing 
methods.  The  same  number  of  independent  tests  were 
performed  with  each  method.  The  mean  strength  is  the 
same  for  both  type  tests,  but  the  scatter  is  larger  for 
test  B.  The  safety  factors  of  an  embankment  constructed 
on  this  foundation  are  identical  for  both  strengtii  deter- 
minations iiiised  on  ttic  average  Su.  However,  the  com[Juted 
pt obab 1 1 1 ty  of  failure  for  method  A is  lower  due  to  the 
smaller  scatter  about  the  mean.  This  case  illustrates 
the  importance  of  random  testing  error. 

The  second  case  demonstrates  the  determination  of 
strength  for  different  subsoils  using  a "perfect"  test 
for  each.  Again,  the  mean  strength  is  the  same  for  both 
s.ubsoi  Is  but  the  scatter  is  different.  The  embankments 
on  both  subsoils  will  have  the  sam.c  safety  factors,  but 
ditli  r<  nl  comfiutod  failure  probabilities.  This  ease 
firei'.ents  the  i mfi  1 i cat  i ons  of  inlu’rent  soil  var  i <ib  i 1 i ty  . 

I i pire  1-3  presents  t;he  importance  of  the  sc.-at.ti'r 
III  eojieetion  fiictor  corre  1 at.  i ons  to  compensate'  tor  bias. 
Till'  sources  of  data  for  correction  factor  A were  I'nibank- 
ment  1 allures  on  soft  clays.  These  opibankment  failures 
were  not  restricted  to  a particular  geologic  deposit. 
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The  sources  of  data  for  correction  factor  B were  retricted 
to  embankment  failures  on  a particular  geologic  deposit 
(eg  Boston  Blue  Clay) . The  mean  correction  factors  based 
on  the  same  number  of  data  are  identical,  but  the  scatter 
in  developing  correction  factor  A was  larger  than  for  B. 

Tf  an  embankment  were  constructed  on  Boston  Blue  Clay, 
tlu’  safely  factor  would  be  identical  regardless  of  which 
corti'ct  ion  factor  wore  applied  to  the  measured  strength. 
llowiv<-r,  applyiiK)  correction  factor  A would  result  in  a 
higher  compuic'd  probability  of  failure  than  correction 
factor  B.  The'  difference  is  due  to  the  uncertainty  in 
the  bias  correction. 

1.4.2  Input  Value  of  Strength.  Once  all  available 
strength  data  are  gathered,  the  engineer  is  faced  with 
the  deicision  of  what  value  of  strength  to  use  in  perform- 
ing his  stability  analyses.  Some  engineers  may  select  to 
use  the  true  mean  of  the  data,  assuming  that  the  safety 
f .ictor  will  .illow  for  tlie  sc.itter  of  the  data.  Sowers 
(l'»70)  recommends  selecting  a value  of  strength  for  which 
2')  fierci.nt  ot  the  data  are  lower  and  75  percent  higher. 
Other  ('ii<}ineers  may  select  the  lowest  observed  strength 
for  design. 

Clf'arly,  only  f»y  start  in<j  with  the  best  e.stimate 
or  mean  strength  (assuming  no  bias  or  <if ter  correcting  it) 
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can  one  arrive  at  the  mean  safety  factor.  If  lower 
strencTth  values  are  assumed  for  design,  the  true  mean 
safety  factor  of  the  embankment  increases.  For  example, 
if  one  designs  an  embankment  for  a safety  factor  equalled 
to  1.5  and  strengths  selected  wore  based  on  Sowers  (1970) 
recommendation,  the  true  mean  safety  factor  could  be 
closer  to  2. 

Because  there  is  no  universally  accepted  method  of 
selecting  the  input  strength,  the  mean  safety  factor  of 
previously  constructed  embankments  is  uncertain.  Also,  a 
set  of  foundation  strength  data  can  have  at  least  three 
different  values  of  input  strength  - the  mean  value. 
Sowers  value  and  the  lowest  value  of  strength.  For  a 
specified  design  safety  factor  and  the  three  different 
values  of  input  strength,  throe  different  embankment 
con  t i cjurat  ions  can  be  designed.  All  three  would  have  the 
same  safety  factor,  but  the  calculated  failure  probabil- 
ities would  bo  quite  different. 

1 . 5 IJncor  ta  i n_ti^  i_r^  the  Load 

1.5.1  Density  and  Configuration  Measurements.  The  soil 
property  needed  for  load  calculations  is  the  total  unit 
wf'icjht  of  the  soil.  This  property  is  potentially  subject 
fo  the  s.mie  sources  of  uncertainty  as  is  strength  - 
namely  bias,  random  testing  error,  and  inherent  soil 
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var  1 .il)  1 1 1 ty  . However,  unit  weiqht  determinations  are 
considered  to  be  unbiased  tests. 

Also  required  for  calculating  the  load  caused  by 
the  unbalanced  mass  of  an  earth  embankment  is  the  config- 
uration of  the  embankment.  For  most  man-made  earth 
structures  this  configuration  is  well  defined  and  varia- 
tions in  configurations  are  small. 

1.5.2.  External  Loads.  The  uncertainty  of  external 
sources  of  loading  also  affects  the  failure  probability. 
These  external  sources  can  be  classified  as  dynamic  and 
static.  Dynamic  loads  include  those  due  to  earthquakes, 
(’xplosives  or  impact.  Both  the  magnitude  and  the 
uncertainty  of  these  loads  can  bo  large.  These  type- 
loads  have  uncertainty  both  in  magnitude  and  in  time  of 
occurance . 

Although  they  may  be  of  ma^or  consequence  for  certain 
slopes,  the  uncertainty  of  loading  is  usually  much  smaller 
than  that  of  resistance.  For  the  general  case  developed 
in  this  report,  tlie  uncertainty  of  the  load  is  necjlected 
in  the  probability  model,  i.e.,  the  load  is  chosen  to  be 
del  I 1 Ml  1 ri  I s t 1 e . 

I . Inllui-nee  of  _Fml)ankJTHTit  Iierigl]2-  Ttie  current  di'sicjn 
appro. icti  doe;:  .illow  for  li-ngth  of  tJie  emliankmcnt  in  one 
respect.  If  noticeably  diffc-rcuit  soil  properties  or 
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foundation  geometry  are  uncovered  along  the  length  of  a 
proposed  embankment,  the  embankment  is  divided  into 
sections  having  similar  subsurface  conditions.  In  fact, 
such  subdivisions  should  also  be  made  when  applying  a 
probability  analysis  to  the  embankment. 

This  subdivision  of  embankment  length  does  not  fully 
tre.it  the  uncertainty  due  to  embankment  length.  Consider, 
for  exampli',  two  embankments  designed  for  a safety  factor 
equ.illcd  to  1.50.  One  embankment  is  one  mile  long,  the 
other  is  ten  miles  long.  Assume  that  the  subsurface 
conditions  were  similar  for  both  embankments  and  that  the 
embankment  length  was  subdivided  when  necessary. 

If  both  embankments  were  designed  in  a similar 
manner  and  all  other  uncertainties  were  equal,  the  longer 
tmibankment  would  have  a higher  probability  of  a failure 
occur irig  somc-where  along  the  embankment  length.  This  is 
true  even  thougfi  the  safety  factors  are  equal.  The 
effects  of  embankment  length  will  be  developed  in 
Chapter  4 of  this  thesis. 

1.7  Advantages  of  a Probability  Model.  There  are  two  main 
re.isons  why  .i  probability  model  should  supplement  the 
current  design  approach.  The  first  reason  is  to  develop 
1 mi  ■ t liodo  1 o<iy  th.it  can  handle  the  limitiitions  of  the 
convent  lon.il  salt'ty  factor  method  of  analysis.  The  second 
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and  most  important  is  to  havo  an  expression  for  an  earth 
embiinkment  failure  probability  that  can  be  combined  with 
the  costs  associated  with  the  structure  to  optimize  design. 

The  previous  section  stated  the  main  limitations  of 
the  current  design  approach.  It  is  the  ultimate  goal  of 
the  probability  model  to  give  a meaning  to  all  of  the 
important  factors  that  influence  the  probability  of 
failure  of  a structure.  The  probability  model  will  not 
minimize  the  role  of  "engineering  judgement".  Rather,  it 
forces  the  engineer  to  further  use  "engineering  judgement." 
in  evaluating  parameters  such  as  soil  variability,  bias, 
and  random  error  which  are  important  for  the  ultimate 
safety  of  a structure. 

Once  the  safety  of  a structure  is  expressed  in  terms 
of  its  probability  of  failure,  a design  optimization 
procedure  can  be  developed.  This  procedure  gives  the 
engineer  a method  of  detc'rmining  what  probability  of 
failure  to  design  for.  In  t tie  past,  design  safety  factors 
havf'  been  set  based  on  intuitive  judc|ement  and  past 
experience.  In  some  case's  t tie  coriseejuences  of  failure 
tiave  boi'n  overlooke'd.  By  comliining  the  consequf'nces  of 
1. 11  lure-  with  ttie  p roli.ib  1 1 i t y of  Ihiilure  and  ttie  construc'- 
tion  ('ost  s a logical  decision  can  be  made  as  to  the 
reejuire'f]  .safe'ty  of  a structure. 
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Consider  the  case  of  an  earth  dam  analyzed  for  the 
end  of  construction  stability.  By  evaluating  the 
possibility  of  f.ii  lures  for  different  assumed  geometrys 
and  by  estimating  the  construction  costs  for  each  design, 
a plot  of  construction  costs  (cc)  versus  probability  of 
failure  (pf)  can  be  made,  see  figurel*4.  The  next  step  is 
to  estimate  the  monetary  consequence  of  a failure  (CF) . 
Then,  {>lot  the  product  of  (CF)  times  (pf)  versus  (pf). 

The  total  I'xpected  costs  (TKC)  is  equal  to  the  sum  of 
(cc)  plus  (CF)  times  (pf).  The  recommended  design 
probability  of  failure  is  that  which  minimizes  the  total 
expected  costs. 

The  above  procedure  assumes  that  the  conseciuences 
of  failure  are  only  of  an  economic  nature.  This  means 
that  the  risks  of  death  or  bodily  injury  and  the  major 
environmental  consequences  s)iould  also  bo  assigned  costs. 

It  can  be  seen  how  the  current  design  method  can 
overlook  tiie  process  of  Ofit  i m i za  t i pn  by  considering  an 
('arth  dam.  An  e.arth  dam  is.  usu.iNy  desi<jned  for  a safety 
factoi'  of  1.00  for  both  the  alter  construction  case  and 
the  steady  seepage  case.  Simil.ir  s<ifety  factors  imply 
similar  failure  probabilities.  However,  the  conse<]uences 
of  a dam  failing  during  steady  seepage  with  a reservoir 
full  of  water  ari'  orders  of  magnitude  larger  than  those 
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of  a dam  fculing  during  construction. 


By  proceeding  with  a design  optimization  approach 
based  on  economic  considerations,  a logical  decision 
can  be  made  as  to  how  safe  a structure  should  be.  The 
probability  model  assumes  that  all  earth  embankments  have 
a calculable  probability  of  failure.  By  recognizing  and 
evaluating  this  probability,  the  engineer  is  able  to  treat 
it  in  a logical  manner. 


' 


i 
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Chapter  2 The  F’robability  Model 

2.1  Genera  1 . As  seen  in  Chapter  1,  the  conventional 
method  of  designing  earth  embankments  has  two  main  limita- 
tions. It  is  the  intent  of  this  chapter  to  provide  an 
answer  to  one  of  them,  i.e.,  its  inability  to  incorporate 
the  degree  of  unreliability  inherent  in  the  decision-making 
required  for  stability  analyses.  The  probability  model 
will  allow  one  to  evaluate  this  unreliability  and  to  in- 
clude it  as  an  input  parameter. 

In  this  chapter  a two-dimensional  model  for  slope 
stability  problems  on  a homogeneous  foundation  is  developed. 
It  is  followed  by  a discussion  on  how  to  evaluate  the  un- 
certainty quantitatively.  Finally,  the  model  is  extended 
to  handle  stability  analyses  with  layered  foundations. 

2.2  Statistical  Tools.  Before  presenting  the  probability 
model,  some  basic  concepts  and  tools  are  needed.  First, 
it  will  be  assumed  that  the  soil  properties  used  in  the 
model  are  normally  distributed.  A normal  distribution  is 
a hell-shaped  curve  giving  the  distribution  of  the  proba- 
bility associated  with  the  different  values  of  the  data. 

The  normal  distribution  of  soil  properties  of  interest  for 
slope  stability  problems  has  been  suggested  by  Lumb  (196b). 

Any  random  quantity  x can  be  represented  by  two 
parameters  — its  mean  value  (X)  and  its  standard  deviation 
(o^)  or  coefficient  of  variation  (V^) . Assume  that  a set 
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of  observations  of  x are  denoted  by  x,  x^  x.  ...,x  The 

1,2,1,  n. 

mean  value  of  X can  be  estimated  by  the  numerical  average 

( * ) 

of  the  n data  points,  i.e., 


X = Xx^/n 


:2. 2-1) 


The  scatter  of  values  around  the  mean  value  is  measured  by 


the  variance  (ox  ) , which  is  estimated  in  the  following 
(*) 

2 

ox  = X — (2,2-2) 


way : 


By  definition,  the  standard  deviation  (o^)  is  equalled  to 
the  square  root  of  the  variance.  Approximately  seventy 
percent  of  all  data  is  contained  within  the  one  sigma 
bounds  (+a  ) . 

- X 

The  coefficient  of  variation  (V  ) of  x . is  obtained 

X 1 

by  normalizing  the  standard  deviation  by  dividing  it  by 


the  mean,  given  as: 

^ o 

= qr  (2.2-3) 

^ X 

The  coefficient  of  variation  and  the  standard  deviation 
would  both  equal  zero  for  data  having  no  scatter  and  would 
increase  with  increasing  scatter  about  the  mean.  Examples 
of  the  above  definitions  are  shown  in  Figure  2-1. 


♦Representing  the;  estimated  variance  of  X with  the  notation 
sx^  has  been  proposed  by  Benjamin  and  Connell  (1970). 

Thc'ir  proposed  notation  for  "true"  or  theoretical  variance 
(o  ^)  was  used  in  this  r;  ^ ort  for  estimated  valiu'S.  This 
notation  was  selected  to  avoid  notational  difficulties, 
since  soil  strength  is  also  expressed  by  tlie  lett('r  S. 
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If  paired  data  are  available  the  correlation 
between  them  can  be  measured  by  calculating  the  co-variance. 


COV(x,y),  estimated  as  follows: 

I ^ 

COV(x,v)  = - L (x.-  X)(y.-  Y) 
n 1 •'  1 

Normalizing  the  co-variance  by  dividing  it  by 

of  the  standard  deviations  of  the  two  data  (X 

the  correlation  coefficient,  (p  ),  given  as: 
xy  ^ 

^ COV(x,y) 


(2.2-4) 

the  product 
and  Y)  yields 

(2.2-5) 


I 


The  correlation  coefficient  can  vary  between  +1  and 
-1.  Both  extreme  values  indicate  perfect  correlation, 
positive  and  negative  respectively,  between  the  sets  of 
data  (Benjamin  and  Cornell,  1970).  That  is,  the  data  will 
plot  as  a straight  line  on  a natural  scale  plot.  A corre- 
lation coefficient  equalled  to  one  means  that  high  values 
of  X imply  high  values  of  y,  and  vice  versa.  A correla- 
tion coefficient  equalled  to  negative  one  means  high  values 
of  X imply  low  values  of  y,  and  vice  versa.  A correlation 
coefficient  equalled  to  zero  indicates  the  data  are  uncor- 
related linearly.  There  might,  however,  exist  some  higher 
order  corrc!lation  between  the  data  that  is  undetected  by 
the  correlation  coefficient. 

The  correlation  distance  of  a property  (6s)  is  a 
measure  of  how  rapid  the  property  changes  in  space.  This 

requires  calculating  the  coefficient  of  correlation  as  a 
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function  of  distance  apart.  Plottiny  the  coefficient  of 
correlation  versus  distance  apart,  the  correlation  dis- 
tance is  defined  as  that  distance  having  a correlation 
coefficient  equalled  to  e ^ (Vanmarcke,  1974).  Examples 
of  these  definitions  are  illustrated  in  Figure  2“2. 

2 . 3 Basic  Probability  Model 

The  probability  model  presented  in  this  thesis  was 
developed  for  the  circular  arc  analysis  of  slope  stability. 
Assuming  a circular  failure  mechanism,  the  safety  of  an 
earth  slope  against  a shear  failure  can  bo  defined  by  a 
safety  factor: 

Resisting  Moment  M„ 

FS  = r ■ iT — — .=  R = sLr  (2.3-1 

Overturning  Moment  — -- — 

M Wa 

o 


where : 

s = mean  shear  strength  along  the  assumed 
failure  arc. 

L = arc  length  of  assumed  failure  arc. 

r = radius  of  curvature  of  the  failure  arc. 

W = weight  of  the  material  within  the  assumed 
fa i lure . 

a = horizontal  distance  from  the  center  of 
rotation  to  the  center  of  gravity  of  the 
failure  mass. 

For  a description  of  this  exfirt'ssion  of  safety  factor, 
se<'  Figure  2-3. 
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The  safety  of  an  earth  slope  against  failure  can  also 

be  defined  by  its  safety  margin  (M)  equalled  to  the  excess 

resisting  moment.  Considering  and  M to  bo  random 

R o 

variables,  the  safety  margin  is  defined  as; 


M = (M^-  - M ) 

R o 


(2.3-2) 


Fiy  expressing  safety  in  terms  of  net  resistance,  the  problem 
of  handling  the  passive  driving  force  is  eliminated.  The 
safety  margin  will  be  the  same  whether  it  is  considered 
as  a positive  resistance  or  as  a negative  load. 

The  probability  of  failure  (pf)  of  an  earth  slope  is 
equated  to  the  probability  that  the  safety  margin  is  less 
than  or  equalled  to  zero,  expressed  as: 

pf  = p (M  < 0)  (2.3-3) 

By  subtracting  the  mean  safety  margin  from  both  sides  of 
this  equality  and  dividing  both  sides  by  the  standard 
deviation  of  the  safety  margin,  the  probability  of  failure 
becomes : 


pf  = p [ 


M - M 


m 


M 

o 


m 


or : 


where : 


pf  = p [ u < -B I 

U = standardized  Safety  Margin  = 


(2. 3-4) 

(2.  3-5) 

M - M 
0 

m 


B = Reliability  Index  = 


M 


m 
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Expanding  Equation  (2.3-5),  the  reliability  index  becomes: 


- M 

= R o 

2 2 h 

(oM^^  + oM^  ) ' 
K O 


(2.  3-6a) 


M„/M  - 1 
R'^  o 


„ aM  T .... 
(~)^  (_°)2  + ( — °)2 


M 


M 


nM 
c 


(2.  3-6b) 


F.S  - 1 

6 = ZZT5 5 (2. 3 -6c) 

(FS  VM^  +VM^  ) ^ 


where  FS  = the  mean  safety  factor 


= the  coefficient  of  variation  of 
the  resisting  moment. 


VM  = the  coefficient  of  variation  of 
o 

the  overturning  moment. 


Assuming  the  resisting  moment  and  the  overturning 
moment  normally  distributed,  the  safety  margin  is  normally 
distributed.  B,  then,  is  a measure  of  the  number  of 
standard  deviations  between  the  mean  safety  margin  and  the 
safety  margin  equal  to  zero  (t.iilure).  By  evaluating  I' , 
the  probability  of  failure  can  be  ohtained  from  a Gaussian 
or  normal  distribution  table.  A <irap)hical  representation 
of  the  reliability  indf?x  and  a pilot  of  pf  versus  B are 
shown  in  Figure  2-4. 
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The  assumption  of  normality  of  the  distribution  of 
the  standardized  safety  marain  can  be  critical.  The  larger 
the  value  of  the  more  critical  the  assumption  of  the 
distribution  type  becomes.  This  is  due  to  the  greater  im- 
portance of  the  tails  of  the  distributions  at  greater  values 
of  p.  Ang  (1972)  concluded  that  the  distribution  sensi- 
tivity is  not  too  important  for  failure  probabilities 
greater  than  one  in  one  thousand.  The  assumption  of  normal 
distribution  can  lead  to  serious  error  for  failure  proba- 
bilities greater  than  10 

Inspecting  equation  (2.3-6c),  it  can  bo  seen  that  the 
probability  of  failure  of  an  earth  embankment  is  not  only 
dependent  on  the  magnitude  of  the  mean  overturning  and 
resisting  moments,  but  also  on  the  uncertainty  involved  in 
obtaining  them.  Figure  2-5  presents  how  the  reliability 
index,  and  hence  the  probability  of  failure  is  affected  by 
the  magnitudes  of  the  coefficients  of  variations  of  the 
resisting  and  of  the  overturning  moments. 

Figure  2-5a  shows  the  relationship  between  the  relii- 
bility  index  and  the  mean  safety  factor  for  different  coef- 
ficients of  variation  of  the  overturning  moment  with  tlu' 
coc'fficient  of  variation  of  the  resisting  moment  etpj.illed 
to  zero.  Since  a .••.traight  line  relationships  exists  betwe('n 
the  two  variables,  an  embankment  can  be  built  to  take  any 
unreliability  in  the  overturning  moment.  This  is  due  to 


the  fact  that  the  magnitude  of  resistance  is  certain.  The 
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graph  also  indicates  that  variability  in  the  overturning 
moment  can  imply  large  variationsin  the  reliability  index 
for  the  same  value  of  safety  factor. 

Figure  2-5b  shows  the  converse  plot  with  VM^  equalled 
to  zero  and  VM  a variable.  This  plot  indicates  that  the 
value  of  the  reliability  index  approaches  a maximum  value 
equalled  to  1/VM  as  a limit.  This  moans  that  for  a given 

A 

VM  It  may  be  impossible  to  design  for  a specified  proba- 
bility  of  failure  by  only  increasing  the  central  safety 
factor.  The  plot  also  shows  that  variability  in  the  resist- 
ing moment  significantly  affects  the  reliability  index. 

The  above  assumption  of  VM  as  constant  for  increasing 

central  safety  factors  is  conservative.  With  an  increase 

in  associated  with  an  increase  in  the  length  of  the 

critical  failure  arc,  VM  will  decrease.  The  maximum  value 

K 

of  the  reliability  .index  equalled  to  1/VM  is  a lower  bound. 

R 

Since  the  uncertainty  in  resistance  is  usually  much 
larger  than  the  uncertainty  in  load,  the  assumption  of  VM^ 
equalled  to  zero  is  often  made.  For  this  case,  the  mini- 
mu.m  allowable  value  of  the  central  safety  factor  can  be 
determined  from  Figure  2-5b  if  the  failure  probability 
(and  hence  f!)  is  specified  and  VM^^  is  fixt'd. 

It  is  interesting  to  note  that  when  the  ce.itral 
safety  factor  is  equal  to  unity,  the  reliability  index  is 
equal  to  zero.  This  value  of  B implies  a probability  of 


failure  of  onChalf  in  the  case  of  the  Gaussian  probability 
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distribution.  A true  safety  factor  of  one  does  not  mean 
that  an  earth  slope  will  fail,  rather  that  the  slope  has  a 
fifty  - fifty  chance  of  survival. 

2 . 4 Load  and  Resistance  Uncertainty 

As  has  been  previously  shown,  the  uncertainties  in 
determining  the  resisting  moment  and  the  overturning  moment 
are  intimately  linked  to  the  probability  of  failure  of  an 
earth  slope.  For  the  general  case,  the  uncertainty  of  the 
overturning  moment  is  due  to  the  scatter  in  the  total  unit 
weight  of  the  soil  and  to  variations  in  embankment  configu- 
ration. For  this  case,  the  variability  of  the  overturning 
moment  is  usually  much  smaller  than  that  of  the  resisting 
moment  and  will  be  considered  equalled  to  zero.  The  equa- 

I 

tion  for  the  reliability  index  now  becomes:  j 

6 = (2.4-1) 

FS  VM„ 

H 

This  section  will  present  a method  for  determining  the 
coefficient  of  variation  of  the  resisting  moment  that 

should  be  used  with  equation  (2.4-1).  , 

Two  different  concepts  of  mean  strength  and  variance 
of  strength  have  to  be  distinguished  from  one  another.  They 
are  the  point  mean  and  point  variance  and  the  spatial  mean 
and  spatial  variance.  The  difference  between  the  two  is 
best  illustrated  with  an  example.  Consider  a "homogeneous" 
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soil  deposit  on  which  twenty- five  strength  measurements 
were  made.  The  moan  and  variance  of  these  data  can  be  cal- 
culated by  equations  (2.2-1)  and  (2.2-2).  Assume  that  the 
mean  and  variance  of  the  data  were  calculated  to  be  500  psf 
and  10,000  psf,  respectively.  The  standard  deviation  equals 
100  psf. 


These  numbers  are  the  best  estimate  and  uncertainty 
of  the  outcome  if  an  additional  test  were  performed.  That 
is,  if  one  additional  test  were  performed,  the  best  estimate 
of  the  resulting  strength  would  be  500  psf.  Also,  there 
would  be  a seventy  percent  chance  that  the  resulting 
strength  would  be  between  400  and  600  psf. 

The  spatial  mean  and  spatial  variance  reflect  the 
best  estimate,  and  the  uncertainty  in  the  best  estimate,  of 
the  average  strength  of  the  soil  deposit.  If  the  strength 
tests  are  uncorrelated,  the  spatial  variance  is: 


os  _ 10 , OOOpsf  _ 


n 2 5 

and  the  spatial  standard  deviation 


4 00  psf 
becomes : 


os  loo  psf  , 

o<s)  = — = tr  = 2 0 psf 

y~n  /~2T' 

The  spatial  mean  is  equal  to  the  point  mean.  Following  up 
on  the  oxamole,  the  best  estimate  of  the  average  strencjth 
of  the  soil  deposit  is  500  psf.  There  is  seventy  percent 
confidence  that  the  strength  of  the  deposit  is  between  480 
and  520  psf. 
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The  measured  strength  at  a point  in  a "homogeneous 


soil  deposit  can  be  expressed  as  a sum  of  contributions: 


S,  \ ^ ^ U 

(meas)  (true)  r b 


(2.4-1) 


where : 
S 


(meas) 


= the  measured  strength  at  a point. 


S,  , = the  true  in-situ  strnegth  at  a point, 
(true) 

^ ^ = the  random  testing  error  at  a point. 

^ = the  bias  or  systematic  error  at  a point, 

The  measured  strength  at  the  i^^  point  is: 

(i)  (i) 


S = S . +€  +£, 

(meas)  (true)  r b 


;2.4-2) 


The  random  error  term  has  a zero  mean  and  a variance,  o €^. 

The  bias  error  term  has  a mean  equal  to£j^  and  a variance 
2 

o€|^.  The  point  average  or  best  estimate  of  strength  at  a 
point  can  be  expressed  as  the  average  of  the  contributions 
to  the  measured  strength: 


g = g + c = g+£ 

(meas)  (true)  ^b  ^b 


(2.4-3) 


The  point  variance  of  the  measured  strength  is  equal 
to  the  sum  of  the  variances  of  its  contributions,  given  as: 


VAR  [S,  1 = + oc  ^ + 0£,^ 

(meas)  ‘ s '-r  b 

(i) 


(2.4-4: 


where:  VAR  [S,  '"'.I  = the  point  variance  of  the 

‘ (meas) ' ^ 

measured  strength  based  on  i independent  tests. 
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I 


2 

= tho  variance  of  the  true  strength,  or  scatter 
due  to  tho  inherent  variability  of  soil. 

The  point  variance  of  the  measured  strength  reflects  the 
confidence  in  the  value  of  strength  obtained  from  a single 
independent  test. 

As  the  number  of  independent  tests  increases,  the 
uncertainty  in  the  effect  of  random  testing  errors  decreases. 
Since  random  testing  errors  fluctuate  about  the  true  moan 
strength,  these  errors  will  be  self-compensating  for  large 
numbers  of  tests.  This  is  not  so  for  the  bias  error.  Test 
bias  is  present  when  a particular  type  of  testing  device 
consistently  measures  strength  either  too  high  or  too  low. 
This  source  of  error  cannot  be  reduced  by  increasing  tho 
number  of  independent  tests,  but  will  be  a constant  source 
of  error  in  trying  to  measure  in-situ  strength. 

The  point  variance  is  not  the  variance  of  interest  for 
probability  models.  What  are  needed  are  the  mc>an  and  t he 
variance  of  the  average  of  all  points  within  the-  region  of 
interest,  or  tht?  spatial  mean  and  the  spatial  variancf  of  r; 
tests  performed.  This  will  give  a measure  of  confidence  1 1' 
the  value  of  averaeje  strength  obtained  by  testing. 

The  spatial  mean  or  the  moan  of  the  spatial  average 


of  all  points  within  the  volume  of  interest  is  equal  to  the 
po i nt  mean , i . e . , 

''S  meas  = • S>true  4 C.  = S 4-  f (2.4-5) 

b b 


where : 


f 


<S>meas 


volume 


S , . dx  dy  dz 

(x,y,z) 


The  proper  evaluation  of  the  spatial  variance,  and 
hence  the  spatial  coefficients  of  variation,  is  a critical 
issue  in  determining  meangful  values  of  failure  probabi- 
lities from  probabilistic  models  of  slope  stability. 

The  spatial  variance  is  shown  in  section  (2.6-4) 

equal  to:  2 2 

- a a 

VAR  [<S>meas)  = o,  + — ^ ^ (2.4-6) 

b n N 

e 

whore:  n = the  number  of  independent  strength  tests 

performed . 

= the  equivalent  number  of  independent  soil 
elements  within  the  dssumed  failure  surface. 
This  number  is  a function  of  the  correlation 
coefficient  of  the  soil  and  the  area  of  the 
failure  surface. 

The  evaluation  and  implications  of  this  spatial  variance  is 
discussed  in  section  (2.6). 

2.5  Model  Uncertainty.  As  pointed  out  in  Chapter  One, 
there  are  two  main  sources  of  uncertainty  due  to  the  select- 
ed method  of  analysis.  The  first  is  the  uncertainty  of  the 
model  (circular  arc  versus  wedge)  and  the  definition  of 
safety  factor.  The  second  is  the  uncertainty  due  to  sim- 
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plifying  assumptions  made  both  for  the  model  and  for  ttio 
problem.  This  model  uncertainty  will  be  handled  in  these 


ways : 

1.  The  probability  model  developed  in  this  thesis 
IS  only  considered  valid  for  circular  arc  type 
failures.  Wedge  methods  of  analysis  are  excluded. 

2.  The  safety  of  an  embankment  is  defined  in  terms 
of  safety  margin  instead  of  safety  factor.  This 
eliminates  the  problem  of  the  location  of  the 
passive  moment  in  safety  factor  definitions. 

3.  The  bias  correction  factors  used  in  this  thesis 
are  developed  from  case  histories  of  embankment 
failures.  The  bias  factor  for  strength  will  also 
include  bias  (and  uncertainty)  due  to  the  method 
of  analysis. 


The  corrected  strength  need  not  be  the  true  in-situ 
strength.  It  is,  however,  the  value  of  strength  that  makes 
the  circular  arc  method  of  analysis  work  — that  is,  pre- 
dict a factor  of  safety  equal  to  unity  for  a failed  embank- 
ment. 


2 . 6 Measurement  of  Resistance  Uncertainty 

Measuring  the  variability  of  the  resisting  moment 

requirc’s  evaluating  the  spatial  uncertainty  duo  to  bias 

corrections,  random  testing  errors,  and  soil  variability. 

These  sources  of  uncertainty  will  be  discussed  separately. 

-51- 


Uncertainty  Due  to  Bias 
The  determination  of  the  means  and 
coefficients  of  variation  of  bias  correction  factors  asso- 
ciated with  different  methods  of  testing  (i.e.,  field  vane 
and  unconfined  compression,  among  others)  is  difficult  to 
determine.  Although  embankments  are  designed  for  a speci- 
fic safety  factor,  the  actual  safety  factor  of  a satisfac- 
torally  performing  embankment  is  unknown.  It  is  only  when 
an  embankment  fails  that  the  true  safety  factor  is  known. 

For  cases  of  embankment  failures,  the  true  shear 
strength  can  be  back  calculated  by  assuming  the  safety 
factor  equal  to  unity.  If  soil  tests  are  available  at 
these  sites,  a bias  term  can  be  determined  by  relating  the 
true  strength  to  the  measured  strength.  Such  a bias  term 
will  include  bias  contributions  from  both  testing  inaccura- 
cies and  model  limitations.  When  a sufficient  number  of 
independent  bias  terms  have  been  developed  for  the  same 
type  strength  test,  an  average  or  mean  correction  factor 
can  be  calculated.  By  calculating  the  coefficient  of  varia- 
tion (V^)  about  this  mean  with  equation  (2.2-3),  the  uncer- 
tainty in  applying  the  mean  correction  factor  is  quantified. 
Such  average  correction  factors  have  been  developed  for 
field  vane  tests  and  unconfined  cempression  tests.  These 
will  be  presented  together  with  the  case  studies. 
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In  developing  a mean  correction  factor  for  a particu- 
lar type  of  strength  test,  care  must  be  taken  as  to  the 
source  of  data.  Ideally,  correction  factors  should  be 
obtained  trom  test  embankments  purposely  constructed  to 
failure,  with  prior  strength  tests  available.  If  the  cor- 
rection factors  are  obtained  from  embankments  that  were 
constructed  to  a reasonable  safety  factor  based  on  prior 
strength  information,  but  failed,  an  erroneous  correction 
factor  would  be  developed.  If  the  strength  tost  used  for 
the  design  of  the  embankment  is  calibrated  to  the  strength 
of  the  failed  embankment,  the  correction  factor  will  equal 
the  inverse  of  the  safety  factor,  neglecting  the  strength 
of  the  embankment.  If  the  bias  correction  for  a given 
embankment  were  greater  than  the  inverse  of  the  safety 
factor,  the  embankment  would  not  fail  and  no  data  on  the 
magnitude  of  the  bias  correction  factor  would  be  obtained. 

The.*  uncertainty  due  to  only  bias  cannot  be  decreased 
by  increasing  the  number  of  independent  strength  tests  at 
an  embankment  site.  Each  individual  test  will  be  biased 
and  the  average  of  all  tests  will  also  be  biased.  The 
coefficient  of  variation  of  the  correction  factor  data 
(correction  factors  obtained  from  a number  of  independent 
case  studies  of  failures)  about  its  mean  value  will  be  a 
measure  of  the  uncertainty  in  applying  that  correction 
factor . 
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[f  the  coefficient  of  variation  of  the  bias  correc- 


tion is  large  and  if  a large  number  of  tests  is  performed 
for  a given  embankment  on  a "homogeneous"  foundation,  the 
unreliability  of  the  bias  correction  will  probably  be  the 
dominant  source  of  uncertainty.  for  this  reaosn  it  seems 
essential  that  su  ch  correction  factors  should  be  updated 
whenever  possible.  This  can  be  done  by  performing  strength 
tests  at  sites  of  embankment  failures  after  the  failure  has 
occurred.  Each  observed  failure  will  results  in  another 
point  for  calibrating  the  different  types  of  tests. 

Although  bias  error  cannot  be  eliminated  by  increased 
testing,  it  can  be  reduced  by  restricting  the  applicability 
of  the  bias  correction.  Attempts  should  be  made  to  develop 
local  correction  factors  and  to  evaluate  their  uncertainty. 
For  similar  type  sampling  and  testing  on  the  same  type  soil, 
the  coefficients  of  variation  of  the  correction  factors 
should  decrease.  Such  an  effort  to  calibrate  a test  like 
the  field  vane  at  a local  level  will  result  in  an  increased 
reliability  of  that  test  for  design. 

2.6.2  Uncertainty  Due  to  Random  Testing  Error. 

This  error  term  is  needed  to  allow  for 
the  error  in  predicting  the  spatial  average  strength  caused 
by  limited  soil  testing.  If  all  of  the  strength  tests  per- 
formed are  far  enough  apart  so  that  they  can  be  considered 
uncorrelated  and  if  the  sum  of  random  errors  is  expressed  as: 

-64- 


(2. 6-1) 


E - £ + £ 

R r,  r-  r 

12  n 


with  the  mean  and  variance  given  below: 


Ed  = 

R r 


OE 


R 


= noc^ 


By  defining 


JR 

n 


(2.6-2a) 
(2. 6-2b) 

(2.6-3) 


then  the  spatial  mean  and  the  spatial  variance  of 
random  testing  error  becomes: 


<Cj_>  = Ej^/n  = n£^/n 


<€.  > = ^ 
r ^ r 


and 


o ' £■ 


a <£  > 

r 


no- 


“2  °T-. 

n^  E. 


(2. 6-4a) 
(2. 6-4b) 

(2.6-5a) 

(2. 6-5b) 


By  ta)cing  the  point  variance  due  to  random  testing  error 
equal  to  that  of  the  measured  strength,  equation  (2.5-5b) 


becomes 


()  <€  > = 
r 


s 

n 


(2.6-6) 


The  cooffici('nt  of  variation  of  the  spatial  nuMn 
shear  strength  available  is: 
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(2. 6-7a) 


/ 2 - , *5  2 

v-v  • <1IVL 

— n ^ r 

<f  ^ 


V 


V-  e > = 

r 


(2. 6-7b) 


The  coefficients  of  variation  of  the  spatial  mean 
shear  strength  due  to  random  testing  error  is  proportional 
to  the  point  coefficient  of  variation  of  the  individual 
data  about  the  mean  and  inversely  proportional  to  the  square 
root  of  the  number  of  tests  performed. 

2 . f).  3 Oncer ta Due  to  Inherent  Soil  Variability 

Even  the  most  uniform  soils  have  some  heterogenity 
of  strength.  The  magnitude  and  rate  of  this  natural  varia- 
bility affects  the  calculated  failure  probability  of  an 
earth  emliankment.  What  is  needed  is  a measure  of  the  un- 
certainty of  encountering  a zone  of  weak  material,  suffi- 
ciently large  for  a failure  to  occur  or  to  adversely  affect 

performance.  Such  a measure  is  the  spatial  variance  due 

2 

to  soil  variability  (o  ^s>).  This  spatial  variance  is  a 
function  of  the  point  variance  of  the  true  strength,  the 
corri'lation  distiincc  of  the  soil,  and  the  area  of  the 
assumed  failure  surface. 

It  can  be  shown  that  the  spatial  variabi 1 ity  of 
striTigth  due  to  soil  heferogen  ity  is  related  to  its  point 
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variance  through  "variance  functions".  First,  consider  the 
concept  of  one-dimensional  variation  of  a soil  property, 


e.g.,  the  variation  of  strength  with  depth  in  a borehole. 

The  spatial  variance  can  be  shown  equal  to 

2 2 2 

o<s  ■ = r , a (2.6-8) 

(h)  s 

where : 

2 

o<s>  = the  spatial  variance  of  strength 
due  to  soil  variability 

2 

r , = the  variance  function 

(h) 

h = the  length  of  the  borehole  considered 

2 

a = the  point  variance  of  strength 

s 


If  the  correlation  coefficient  of  strength  between 
two  points  varies  as  a function  of  the  distance  between 
points  in  the  following  way: 


^s^S2  ^(h) 


2.5-9) 

where: 

I'is  = the  correlation  distance  of  the  strength 
of  soil  considered. 

Then,  Mie  variance'  function  will  equal: 

-hV6s^ 

/■ir 


r 


(h) 


iSs  \ h , , h 

h“  I 7s  • ‘7^, 


) . -I 


(2.6-10) 


where' : 


2 X _ p 

'l'(x)  = the  "';rror  Function"  = ,-  fe  ^ ^ du 
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A plot  of  the  variance  function  versus  h/ 6s  are  shown  in 

Figure  2-6.  The  inverse  of  the  variance  function  is  defined 

as  the  equivalent  number  of  independent  soi.l  ('lements  (n  ). 

e 


n = 
e 


(2.6-11) 


[h) 


If  h is  much  larger  than  the  correlation  distance, 

r ^ reduces  to : , •• 

_ /n  6s  1 


(h) 


--f?  .6-12) 


If  h is  much  smaller  than  the  ccif relation  distance,  !' 


r-  P 


(h) 


etjuals  unity.  This  implies  that  the  points  within  the  length 
h are  fully  correlated  and  that  the  spatial  variance  is 
‘equal  to  the  point  variance. 

The  concept  of  two  dimensional  spatial  variance  of 
strength  over  a rectangular  area  can  also  be  related  to  the 
point  variance  through  variance  functions.  The  two-dimen- 
sional variance  function  is  equalled  to  the  product  of  the 
orthogonal  variance  functions  of  the  two  lengths  of  the  rec- 
tangle (h,{,).  The  spatial  variance  can  be  expressed  as: 

(2. 6-13a) 


o ^ < s > = r ^ ( h , ?,  ) o ' 

s 

u='-'s>  = I’Mh)  (£  )o' 


(2. 6-13b) 


where  i'^(h)  and  r^(f  ) are  the  variance  functions  in  the  h 
and  ^ directions,  respectively.  Again,  if  h and  ? arc  much 
larger  than  the  correlation  distances  in  their  respective 

directions,  ('*' s.  and  ^a.),  then: 

h I 
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r'  {h,  9 


8h 


and : 


th 


e TT6S,  iSs, 


(2.6-14) 


(2. 6-15) 


h'-S 

The  concept  of  one  and  two  dimensional  strength  variability 
due  to  soil  heterogenity  is  illustracted  in  Figure  2-7. 

The  equivalent  number  of  independent  soil  elements 
in  a rectangular  area  of  soil  is  proportional  to  the  area 
of  the  rectangle  and  inversely  proportional  to  the  product 
of  the  orthogonal  correlation  distances.  This  relationship 
implies  that  as  the  correlation  of  strength  increases,  the 
lilcelihood  of  encountering  a large,  uniformly  wea)c  zone  of 
soil  also  increases.  But  as  the  area  of  the  rectangle  in- 
creases, the  likelihood  of  such  a weak  zone  existing  over 
the  whole  rectangle  decreases. 

A simplified  approach  to  estimating  the  equivalent 
number  of  independent  soil  elements  in  an  assumed  failure 
plane  (3-dimensional)  begins  with  an  idealization  of  the 
failure  plane.  Such  an  idealization  allows  one  to  deal  with 
the  three  principal  correlation  distances  in  a straight 
forward  manner.  The  failure  plane  is  modeled  as  a series 
of  rectangular  surfaces  as  shown  in  Figure  2-8.  Consider- 
ing the  cross  sectional  view,  the  length  of  the  failure 
surface  and  the  depth  of  the  failure  are  preserved.  The 
length  in  the  third  dimension,  or  length  of  embankment 
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alonq  the  centerline,  is  selected  as  three  tines  the  base 
width,  D.  This  serves  as  a crude  estimate  of  the  third 
dimension  to  satisfy  the  plane  strain  assumption  of  the 
ana lys i s . 

The  surface  areas  in  the  three  pr incipa 1 planes  are 
then  combined  into  three  rectangular  areas.  The  equivalent 
number  of  independent  elements  is  calculated  for  each 
area  by  the  equations  presented  in  the  previous  section. 

The  sum  of  those  calculated  elements  is  an  approximation 
of  the  number  of  independent  soil  elements  for  a three- 
dimensional  failure  plane  . By  determining  the  point 

variance  of  strength,  the  spatial  variance  will  be  equalled 
to : 


2 oi 

^ <S>  = ^ 

N 

e 


(2. C-16a] 


with  a spatial  coefficient  of  variation  equal  to: 


V' 

s 

<s>  — 

e 


(2. 6-16b) 


2.6.4  Combined  Uncertainty 

The  variances  of  bias  correction, random  error,  and 
soil  variability  can  be  added  together  to  give  the  total 
variance  in  the  prediction  of  strength: 

VAUI<3>  1 = cl’  n 0^ 


meas 


€ b>  + f r > + < s> 


(2. 6-17a) 


or : 
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= ' moas>  ’ »£b  * - * ^ (2.6-17b) 

n N 

e 

The  spatial  coefficients  of  variation  of  the  bias 
correction  factor,  random  error,  and  inherent  soil  varia- 
bility can  be  combined  by  taking  the  square  root  of  the 
sum  of  the  squares  of  the  undivided  coefficient  of  varia- 
t ion : 

V-S>  = (Vb^  + (2.6-18) 

n Ng 


where : 


V<S > is  the  spatial  coefficient  of  variation  of 
strength . 

This  is  a measure  of  uncertainty  equal  to  the  uncertainty 
in  the  resisting  moment  for  the  case  of  an  embankment  on  a 
homgeneous  foundation. 


2 . 7 Probability  Model  for  Layered  Foundations 

Extending  the  general  probability  model  to  handle 
layered  soils  introduces  a second  aspect  of  soil  variabil- 
ity into  the  problem.  That  is  the  variability  in  the  thick- 
ness of  the  different  layers  encountered  through  subsurface 
exploration.  The  thickness  of  a layer  of  soil  is  measured 
by  a direct  observation  of  samples  from  boreholes.  This 
type  of  measurement  is  assumed  to  be  unbiased  and  free 
from  random  error.  The  spatial  variance  of  a layer  thick- 
ness is  a function  of  the  point  variance  of  layer 
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thickness,  (o^^),  and  the  correlation  distance  of  layer 
thickness,  and  the  assumed  area  of  an  embankment  failure. 

The  point  variance  can  be  obtained  by  calculatinq 
the  variance  of  observed  thicknesses  from  a series  of 
boreholes.  The  point  variance  will  give  an  indication  as 
to  what  degree  the  layer  thickness  will  vary  along  the 
length  of  an  embankment.  The  spatial  variance  of  layer 
thickness  also  can  be  related  to  the  point  variance 
through  variance  functions: 

0^  = (Z)  rUw)a\  = ^ (2.7-1) 

o 

whore : 

Ne  = 

rM  5.)  r"  (w) 


I?  = the  length  of  the  layer  parallel  to  the  em- 
bankment centerline  included  in  the  assumed 
failure  surface; 

w = th€“  width  of  the  layer  contained  within  the 
assumed  failure  arc. 

r^(f.),r’(w)  = the  horizontal  variance  functions 

parallel  to  and  fierpendicular  to  the  embank- 
ment centerline,  respectively, 

= the  equivalent  number  of  independent  layer 
thicknesses  within  an  assumed  failure  plane. 

These  variance  functions  will  be  assumed  to  be  of 
the  form: 
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the  np.itial  coefficient  of  variation  of  layer  thickness 
becomes : 

V = V,f//N^  (2.7-5) 

Till'  implication  of  .spatial  variability  of  layer 
thickness  are  illustrated  in  Fiqure  2-9.  Both  embankment 
foundation  A and  B contain  a layer  of  weak  soil.  The  mean 
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layer  thickness  and  the  point  coefficient  of  variation  of 
layer  thickness  arc  the  same  for  both  cases.  The  correla- 
tion distances  for  case  A are  much  greater  than  those  for 
case  B.  The  extent  of  the  assumed  failure  surface  is  also 
shown  in  Figure  2-9  . For  case  A,  the  spatial  coefficient 
of  variation  will  be  approximately  equalled  to  the  point 
coefficient  of  variation.  For  case  B,  the  spatial  coeffi- 
cient of  variation  will  approach  zero. 

These  results  are  reasonable.  As  the  correlation 
distances  increase  or  as  the  area  of  the  failure  decreases, 
the  probability  of  enco  intering  a thick,  weak  layer  over 
the  entire  failure  area  increases.  As  the  correlation  dis- 
tances of  layer  thickness  decreases  or  as  the  area  of  the 
failure  increases,  the  chances  of  encountering  such  a thick 
weak  layer  over  the  entier  failure  area  decreases.  If  the 
layer  thickness  fluctuates  rapidly  compared  to  the  failure 
area,  an  averaging  effect  takes  place. 

Once  the  spatial  coefficient  of  variation  of  soil 
strength  and  layer  thichness  are  measured,  the  coefficient 
of  variation  of  the  resisting  moment  can  be  determined. 
Based  on  the  section  in  Figure  2-10,  the  mean  resisting 
moment  of  a failure  arc  passing  through  n layers  of  soil  is: 

M = r EI.s.  (2.7-6) 

K ^ 
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The  variance  of  the  resisting  moment  depends  on  the 
variances  of  and  the  correlations  among  the  products,  f, . s . 
Assuming  that  strength  and  length  are  both  perfectly  cor- 
related within  the  i^^  layer  (which  would  be  reasonable  if 
the  horizontal  correlation  distances  are  relatively  large) 
and  that  they  are  both  uncorrelated  from  layer  to  layer, 
the  variance  of  the  resisting  moment  can  bo  expressed  as: 


? 2 2 — 2—  2 
oMj^  ■ = 7.  1 V's.>  ) (P.  s^  ) 

+ 21(V<t.>^  + V-s.>^)^'  • I.i. 

1 1 11 


(2.7-7) 


V-  J.^>  and  V<s^>  are  the  spatial  coefficients  of  variation 
of  layer  thickness  and  soil  strength. 

and  s^  are  the  mean  thickness  and  shear  strength  of  the 

. th 

1 layer. 

The  coefficient  of  variation  of  the  resisting  moment 
is  obtained  by  dividing  both  sides  of  equation  (2.7-7)  by 
the  mean  resisting  moment  squared,  or: 


7.  (V<t.>^+  V"s.-^)  ^^i  ^i  ^ + 2MV^<t.>  + V^<s.  ■ ) ^‘' 


n 


{YA.s.}^ 

11 


{Z.s.) 

11 

(j;i^s. ) 


(2.7-8) 


The  derivation  of  this  equation  is  presented  in 
Apjjendix  A . 
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Chapter  3 Case  Studies 


3 . 1 General 

The  application  of  the  probability  model  presented 
in  the  previous  chapter  is  illustrated  with  two  examples. 

In  this  chapter,  a prcbabal  istic  analysis  is  made  fer  tv.o 
different  embankment  case  studies.  The  results  obtained 
from  the  orobabal istic  analysis  are  compared  to  those  from 
a conventional  stability  analysis. 

3 . 2 Fore  River  Test  Section 

3.2.1  Convcnticral  Analysis 

The  embankment  analyzed  in  this  chapter  was  a test 
section  constructed  for  the  Maine  State  Highway  Commission. 
The  Fore  River  Test  Section  (FRT)  was  designed  in  order  to 
evaluate  the  efficiency  of  sand  drains  for  the  proposed 
I-20S  hicjhway  {Haley  and  Aldrich,  1967).  A plan  map  for 
this  test  section  is  shown  in  Figure  3-1. 

The  240  foot  square  tost  section  was  constructed  at 
Portland,  Maine,  over  a tidal  mud  flat.  The  natural  ground 
surface  at  the  site  of  the  test  section  was  approximately 
horizontal  and  at  El.  -3  feet,  msl . The  first  phase  of 
construction  required  placing  granular  fill  to  the  sectio. 
shown  on  Figure  3-2. 

Sand  drains  then  were  to  be  installed  through  the 
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first  phase  of  construction  into  the  foundation  clay.  This 
was  to  be  followed  by  additional  granular  fill  to  El.  +20 
ft.  However,  shortly  after  construction  of  the  first  lift, 
an  undrained  shear  failure  occurred  on  the  North  slope  of 
the  test  section. 

The  foundation  for  the  test  section  consisted  of  a 
layer  of  soft  gray  organic  silty  clay  with  shells,  sand 
lenses  and  pockets,  and  wood  chips.  This  clay  layer  ex- 
tended to  El.  -30ft,  where  silty  sand  and  gravel  was  en- 
countered. 

The  underlying  sand  deposit  had  a 5 foot  artesian 
head.  The  top  of  the  sand  was  considered  to  be  the  limit- 
ing depth  of  failure. 

The  organic  clay  was  tested  both  before  and  after 
the  slide  occurred.  Before  the  failure,  a number  of  boring 
wore  made  and  field  vane  tests  were  performed.  Four  undis- 
turbed borings  also  were  made  with  undisturbed  samples  ta- 
ken at  regular  intervals.  LaJioratory  testing  of  those  sam- 
ples included  visual  classification,  moisture  content  de- 
termination, Atterbcrg  limits,  and  unconfined  compression 
strength  tests  (UC) . 

Within  two  weeks  after  the  slide  occurred,  six  addi- 
tional borings  were  made  through  the  failed  embankment. 
Field  vane  tests  were  performed  in  three  of  these  borings, 
tlnconfined  compression  tests  (UC)  and  unconsol  idat  eei , 


undrained  triaxial  tests  (UUC)  were  performed  on  undis- 
turbed samples  recovered  from  the  other  three  borings. 

Since  the  failure  occurred  in  an  essentially  undrained 
state,  these  strength  tests  will  be  considered  as  measures 
of  the  undisturbed  strength  prior  to  construction,  supple- 
menting the  original  strength  data. 

The  field  vane  strength  data  from  tests  performed 
both  before  and  after  the  slide  are  plotted  in  Figure  3-3. 
Similar  plots  for  the  unconfined  compression  test  data  are 
shown  in  Figure  3-4.  None  of  these  plots  indicate  a defin- 
ite change  in  strength  with  depth;  therefore,  it  was  as- 
sumed that  the  strength  in  the  organic  clay  layer  was  con- 
stant with  depth.  For  the  cases  of  field  vane  tests  and  un- 
confined compression  tests,  the  additional  testing  subse- 
quent to  failure  did  not  change  the  mean  strength  signifi- 
cantly . 

Conventional  circular  arc  and  sliding  wedge  slope 
I stability  analyses  were  performed  on  the  test  section,  con- 

‘ structed  to  El  +10.  The  four  mean  strengths  from  Figures 

3-3  and  3-4  were  used  for  compa r i son . The  critical  failure 
arc  for  the  test  section  is  shown  in  Figure  3-2.  The  cal- 
culated saff'ty  factors  baseil  on  field  vane  data  and  uncon- 
f 1 ne-d  compression  data  were  then  corrected  for  testincj  and 
mode  line)  bias.  This  was  done  by  multiplying  them  by  cor- 
rection factors  proposed  for  the  respective  type  of  strength 
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By  back  calculating  the  in  situ  shear  strength  at 
sites  of  embankment  failures  and  by  comparing  it  to  the 
measured  vane  strength,  Bjerrum  developed  a correction 
factor,  u: 


= u 


s.,  (in  situ) 


(3.2-1) 


Su (FV) 


which  is  a factor  of  the  plasticity  index  (PI)  of  the  soil. 
This  correction  factor  has  been  updated  based  on  additional 
case  studies  by  Milligan  (1972),  Ladd(1973),  and  Flaate  and 
Preber ( 1 974 ) . The  in  situ  strength  is  calculated  as  the  pro- 
duct of  the  measured  strength  and  u,  knowing  the  plasticity 
index  of  the  soil. 

The  resisting  moment  for  the  FRT  test  section  is  al- 
most directly  proportional  to  the  shear  strength  of  the 
organic  clay  layer.  Then,  multiplying  the  safety  factor 
based  on  the  measured  strength  by  the  revised  Bjerrum  fac- 
tor yields  the  best  estimate  of  the  true  (FV)  safety  fac- 
tor. The  revised  Bjerrum  correction  factor  is  shown  in  Fi- 
gure 3-5. 


A similar  correction  factor  was  proposed  by  Wu(1974) 
for  soft  to  medium  clays.  By  relating  the  in  situ  shear 
strength  obtained  from  embankment  failures  to  the  shear 
strength  moasurt?d  by  unconfined  compression  tests,  tlv  Wu 
correction  factor  is  defined  as: 
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(3.2-2) 


( in  situ) 

Su(UC) 

Again,  tiie  best  estimate  of  the  true  (UC)  safety  factor 
at  the  P’RT  test  section  is  the  product  of  the  measured  safe- 
ty factor  and  the  Wu  correction  factor.  This  correction 
factor  developed  from  18  case  studies  has  a mean  value  of 
1.10.  The  Wu  correction  factor  is  shown  in  Figure  3-6. 


The  calculated  and  the  corrected  safety  factors  of 
the  FRT  test  section  for  the  four  mean  strengths  arc  shown 
in  Table  3-1. 


CASE 

Su 

FS 

CF  ((3PI  = 35) 

FS  X CF 

FV ( before 

failure) 

520 

2.04 

0.85 

1.73 

FV (before 

& after  failure) 

520 

2. 04 

0.85 

1.73 

UC  ( before 

failure ) 

210 

0.82 

1.  10 

0.90 

UC (before 

& after  failure) 

205 

0. 80 

1.10 

0. 88 

STABILITY  OF  FRT  TEST  EMBANKMENT 
Table  3-1 


The  resulting  safety  factors  indicate  that  the  testinej 
after  failure  measured  the  same  undrained  strength  as  those 
performed  before  the  failure;  therefore,  the  additional 
testing  did  not  change  the  computed  safety  factor.  Also, 
the  correction  factors  brought  the  safety  factors  based  on 
field  vane  and  on  unconfined  compression  data  closer  to- 
gether. However,  the  corrected  results  liased  on  field  vane 
data  indicate  satisfactory  performance  (FS  =73  percent 
higti);  tlie  corrc'ct  <'d  rosultr.  ))ased  on  unconfined  compres- 
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sion  data  still  show  the  safety  factor  about  10  percent 
too  low. 

Before  the  occurrance  of  the  failure,  it  was  diffi- 
cult to  decide  what  type  of  strength  was  more  reliable  for 
this  test  section.  A number  of  sources  of  bias  control  the 
reliability  of  the  measured  strength  from  both  types  of 
tests  (Haley  and  Aldrich,  1967).  For  the  unconfined  com- 
pression test,  sample  disturbance  and  inhomogeniety  due  to 
sand  lenses  and  shells  would  cause  the  measured  strength 
to  he  too  low.  On  the  other  hand,  the  applied  stress  sys- 
tem and  the  rate  of  strain  to  failure  would  cause  the  mea- 
sured strength  to  be  too  high.  The  net  result  would  pro- 
bably be  strengths  a little  too  low.  This  trend  is  re- 
flected by  the  Wu  correction  factor. 

For  the  field  vane  test,  sample  disturbance  would 
result  in  too  low  strength  measurements  - but  not  as  low 
as  for  unconfined  compression  tests.  This  would  be  ex- 
pected, since  the  vane  test  is  performed  in  situ,  minimi- 
xinej  sample  disturbance.  The  effects  of  the  applied  stress 
system  ,ir<.'  generally  unknown.  Both  the  applied  rate  of 
s.train  (lO'i  [>er  second)  and  the  presence  of  shells,  peat 
fiber.s  and  wood  ch  i fjs  would  caus(>  the  measured  strength  to 
bo  too  high.  The  net  rtnuilt  for  the  field  vane  test  would 
probal)ly  be  strengths  too  high.  This  trend  is  reflected  in 
the  r(‘visc.*d  Bjerrum  correction  factor. 
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The  design  engineer  is  faced  with  the  decision  of 
what  safety  factor  to  select.  Such  a decision  is  conven- 
tionally based  on  experience  and  "engineering  judgement". 


The  following  section  will  examine  the  implications  of 
these  safety  factors,  considering  the  uncertainty  involved 
in  their  determination.  > 

3.2.2  Probabalistic  Analysis 

The  corrected  safety  factors  presented  in  Table  3-1 
are  the  mean  or  central  safety  factors  for  the  four  cases 
based  on  the  available  strength  information.  In  order  to 
assess  the  failure  probability  based  on  these  different 
states  of  knowledge,  it  is  necessary  to  evaluate  the  un- 
certainty associated  with  the  number  and  types  of  test  made 
for  each  case.  As  previously  described,  the  sources  of  un- 
certainty are  in  the  bias  correction  factor,  random  test- 
ing error,  and  inherent  soil  variability.  These  sources 
will  be  looked  at  seperately. 

BIAS  Figures  3-5  and  3-6  present  the  revised  Bjerrum  cor- 
rection factor  and  the  Wu  correction  factor,  repsectively . 
The  revised  Bjerrum  correction  factor  is  the  best  fit  or 
mean  of  the  twenty-eight  available  case  studies.  The  assump- 
tion was  made  that  the  standard  deviation  about  this  mean 
If.  .1  constant.  This  implies  that  the  coefficient  of  varia- 
tion ol  the  correction  factor  is  a variable  and  is  a func- 
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tion  of  the  plasticity  index  of  the  soil.  The  standard 
deviation  of  the  revised  Bjerrum  correction  factor  is 
estimated  at  0.15.  This  correction  factor  and  its  esti- 


mated one  sigma  bounds  are  shown  in  Figure  3-5. 

The  organic  clay  foundation  for  the  FRT  test  section 
had  an  average  plasticity  index  equal  to  35  percent.  This 
plasticity  index  implies  that  a mean  correction  factor  of 
0.85  and  a coefficient  of  variation  (Vb)  due  to  bias  of 
0.15/0. .35  or  0.18  for  the  field  vane  test  results. 

The  Wu  correction  factor  for  unconfined  compression 
tests  is  not  considered  to  be  a function  of  soil  property. 
The  mean  correction  factor  is  a numerical  average  of  the 
correction  factors  of  the  eighteen  available  case  studies 
and  is  equal  to  1.10.  The  coefficient  of  variation  (Vb) , 
of  the  Wu  correction  factor  is  equal  to  0.30. 

RANDOM  TESTING  ERROR  The  coefficient  of  variation  due 
to  random  error  was  shown  to  be  proportional  to  the  point 
coefficient  of  variation.  Vs,  of  strength  and  inversely 
proportional  to  the  square  root  of  the  number  of  indepen- 
dent strength  tests.  The  point  coefficient  of  variation  and 
number  of  tests  [jerformed  for  the  four  mean  strengths  used 
in  analyzing  the  FRT  tost  section  art?  shown  in  Table  5-2. 
Also  shown  in  this  table  arf>  the  s[5atial  coefficients  of 
variation  (Vs//n  ),  for  the  four  cases  analyzed. 
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The  resulting  coefficients  of  variation  indicate  that 
the  number  of  tests  performed  prior  to  failure  was  suf- 
ficient for  both  field  vane  and  unconfined  compression 
tests.  Increasing  the  number  of  independent  tost  neither 
changed  the  mean  strength  significantly  nor  reduced  the 
uncertainly  related  to  random  testing  error  measureably. 

The  results  also  show  that  although  the  standard  de- 
viation of  strength  based  on  field  vane  testing  is  more  than 
twice  that  based  on  unconfined  compression  tests,  the 
normalized  standard  deviationsor  the  point  coof f icients oi 
variation  are  approximately  equal.  The  spatial  coefficients 
of  variation,  (Vs/v'^rT” ),  due  to  random  testing  error  are 
twice  as  high  for  the  unconfined  comression  tests  as  for 
the  field  vane  tests.  This  is  a result  of  the  greater  num- 
ber of  field  vane  tests  performed. 

INHERENT  SOIL  VARIABILITY  The  coefficient  of  variation  duo 
to  inherent  soil  variability  was  shown  to  be  proportional 
to  the  point  coefficient  of  variation  of  strength  and  in- 
versely proportional  to  the  square  root  of  the  equivalent 
number  if  independent  soil  elements.  The  point  coefficient 
of  variation  of  strength,  (V^)  is  shown  on  3-2 

for  the  four  moan  strengths  used. 

The  number  of  equivalent  soil  elements  was  obtained 
by  first  approximating  the  surface  area  of  the  failure  sur- 
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face  With  perpendicular  rectangles.  The  critical  failure 
arc  was  simplified  into  a rectangle,  30  feet  high  and  80 
feet  along  the  base.  The  height  was  selected  equal  to  the 
maximum  depth  that  the  critical  failure  arc  passesd  into 
the  organic  clay  layer.  The  base  width  was  selected  equal 
to  the  difference  in  the  failure  arc  length  and  twice  the 
rectangle  height. 

The  axial  length  of  the  critical  failure  surface 
was  set  equal  to  three  times  the  base  width,  or  240  feet. 
This  length  was  selected  to  preserve  the  assumption  of 
"plane  strain".  The  rectangles  in  the  three  principal 
planes  that  comprise  the  simplified  failure  surface  are 
unfolded  and  shown  in  Figure  3-7. 

The  number  of  equivalent  soil  elements  are  a func- 
tion of  the  dimensions  of  these  rectangles  and  also  o. 
the  correlation  distance  in  the  three  principal  directions. 
A detailed  study  of  the  correlation  of  soil  properties 
is  currently  underway  at  M.I.T.  Although  correlation  dis- 
tances of  soil  strength  have  not  been  investigated  to  date, 
thof-.r-  of  related  material  properties  have  been  measured. 

Based  on  results  of  studies  on  total  unit  weight, 
moisture  content,  and  maximum  past  pressure,  correlation 
distances  of  strength  equal  to  5 feet  vertically  and  100 
feet  horizontally  were  selected.  Using  these  typical  val- 
ues as  estimates  of  the  correlation  distances  of  soil 
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strength  at  the  FRT  test  section,  the  number  of  equivalent 
soil  elements  contained  within  the  simplified  failure  sur- 
face is  determined  from  Figure  2-6. 

The  number  of  equivalent  soil  elements  is  2 for  the 
base  (L  X D)  and  11  for  the  four  sides  (2LH  + 2HB) . The 
total  number  of  equivalent  elements  (Ne)  is  13.  This  num- 
ber is  independent  both  of  the  number  of  tests  and  the  type 
of  tests  performed.  Using  an  average  point  coefficient  of 
variation  for  the  four  cases  equal  to  0.21,  the  spatial 
coefficients  of  variation  duo  to  soil  heterogeniety  is 
eijual  to  0.06. 

A summary  of  the  data  pertinent  to  a probabalistic 
study  of  the  FRT  test  section  is  presented  in  Table  3-2. 

An  important  result  is  that  the  uncertainty  of  the  bias 
correction  factor  (Vb) , as  measured  by  the  coefficient 
of  variation  of  the  bias  correction,  is  the  dominant  source 
of  uncertainty.  This  reflects  the  dilemma  of  the  design  en- 
gineer . 

As  a result  of  the  quantity  of  strength  tests  per- 
formed, the  design  engineer  is  not  much  concerned  with  tlic 
question  of  what  is  the  average  measured  strength.  However, 
he  is  more  concerned  with  the  bias  correction  factor,  or 
which  test,  if  either,  is  measuring  the  in  situ  strength 
correctly  after  applying  a mcMn  bias  correction  factor. 

Thi'  total  uncc'rtainty  in  predicting  the  spatial 
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average  shear  strength  for  the  FRT  test  section  can  be  cal- 
culated wi th  the  equation  (2.6-18).  For  the  field  vane  me- 
thod of  testing: 

Vs  = Vj^j.  = /0.15-’  + 0.06^  + 0.02^  = 0.16 

and  for  the  unconfined  compression  testing  method: 

Vs  = V^r-  = /o73  0 ^ + 0.06'’  + oT'm^  = 0.31 
The  fractional  contribution  of  the  bias  term  to  the  total 
variance  is  88  percent  and  94  percent  for  the  field  vane 
test  and  the  unconfined  compression  test,  respectively. 

rJow  that  both  the  mean  safety  factors  and  the  un- 
certainties in  the  resisting  moment  are  known,  the  relia- 
bility index,  3,  can  be  calculated  from  equation  (2.4-1). 
For  the  field  vane  method  of  testing,  3 equals  2.64  regard- 
less of  which  amount  of  tests  is  perform.ed.  This  value  of 
3 suggests  a failure  probability  of  0.4  percent.  For  the 
unconfined  compression  testing  method  and  a mean  safety 
factor  of  0.89,  3 equals  0.40.  This  value  of  3 loads  to 
an  estimated  failure  probability  of  65  percent.  The  gra- 
phical representations  of  the  reliability  indices  and 
failure  probabilities  of  both  methods  of  testing  are  shown 
in  Figure  3-8. 

The  calculated  failure  probaliility  for  the  uncon- 
fined testing  method  is  high  and  would  gem-rally  b(-  con- 
sidered unacceptable.  On  tlie  other  hand,  the  probability 
of  failure  based  on  field  vane  testing  is  relativf'ly  small 
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and  may  be  considered  as  an  acceptable  risk.  The  fact 
that  the  test  section  actually  did  fail  illustrates  that 
even  with  a high  safety  factor,  an  embankment  is  still  sus~ 
ceptable  to  failure. 

Also  shown  in  Figure  3-8  are  plots  of  failure  probabil- 
ity versus  safety  factor,  assuming  that  the  uncertainty 
in  the  resisting  moment  is  constant. 

3 . 3 Atchafalaya  Basin  Tost  Section 

3.3.1  Conventional  Analysis 

In  order  to  further  illustrate  the  probability  mo- 
I del,  the  stability  of  a flood  levee  is  studied.  The  levee 

in  question  is  part  of  the  Atchafalaya  Basin  protection 
system.  A pair  of  guide  levees,  each  approximately  100 
miles  in  length  were  erected  form  Morganza,  Louisiana  to 
the  Gulf  of  Mexico  to  help  channel  flood  waters  of  the 
Mississippi  River.  By  diverting  the  excess  flood  waters 
■ down  the  Atchaflaya  Basin  Floodway,  the  area  surrounding 

the  lower  280  miles  of  the  Mississippi  River  is  protected 
from  flooding. 

The  construction  of  those  levees  began  in  the  1930 's. 
Due  to  stability  and  deformation  problems,  large  section:; 

h ■ 

. 

ot  t ti('  levees  are  pre:;ent  ly  lielow  design  ijradc.  For  this 

1 

j riMSon  and  because  the  design  flowline  had  been  increased 

I to  allow  for  siltation  in  the  Basin,  several  test  embank- 
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rnents  were  constructed.  These  test  sections  straddled  the 
existing  embankments  and  were  built  to  investigate  design 
alternatives  for  future  construction.  Two  1500  ft  long 
test  embankments,  Test  Sections  II  and  III,  were  adjacent 
sections  that  only  differed  in  their  design  safety  factors. 
Test  Section  II  was  designed  for  a safety  factor  equal  to 
l.i;  Tost  Section  III  for  a safety  factor  equal  to  1.3 
(user;,  1968).  Test  Section  III  (TS  III)  is  analyzed  in  this 
repor  t . 

Before  construction  of  the  e.xisting  levee  began,  the 
ontinal  gournd  surface  was  approximately  horizontal  and 
at  El.  +imsl.  The  foundation  materials  were  very  soft  to 
medium  clays  to  El.  -120msl.  This  thick  clay  layer  was 
underlain  with  sand.  Geologically,  the  clay  is  divided 
into  three  main  deposits.  From  ground  surface  to  El.  -25msl 
is  a poorly  drained  swamp  deposit.  Very  high  moisture  con- 
tent is  encountered  between  El.  -5  and  -25ft.  Below  El. 
-25ft,  the  clay  consists  of  lake  deposits  with  a layer 
ot  well  drained  swamp  deposit  between  El.  -75  and  -90ft. 

The  existing  levof'  was  initially  constructed  with 
mat<,'rial  excavated  from  an  ada  jeent  borrow  [iit.  The;  pit 
was  locati’d  af^prox  ima  tod  1 y 300ft  to  the  floodway  side  of 
th<‘  levee  centerline.  Additional  lifts  ensued  until  l;y 
1965  the  levee  was  at  El  +17ft  and  th(^  borrow  pit  at  El. 
-20ft.  Approximately  6 ft  of  settlement  under  the  exis- 
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tinq  levco  centerline  has  been  observed  from  the  center- 
line  soil  borings  made  in  1964. 


H 

I 


The  sequence  of  construction  of  TS  III  is  shown  in 
Figure  3-9.  Enlarging  the  existing  levee  involved  five 
basic  steps: 

1.  Construct  cast  fill  dikes  on  both  sides  of  the 
existing  levee  to  serve  as  retaining  dikes  for 
the  hydraulically  placed  berms. 

2.  Place  semi-compacted,  hauled  fill  within  the 
main  levee  section  to  El.  +17msl. 

J.  Pump  in  the  hydraulic  fill  berms. 

^ 4.  Place  semi-compacted , hauled  fill  to  the  interim 

gross  grade. 

5.  At  a later  date,  place  semi-compacted,  hauled 
fill  to  the  final  gross  grade. 

The  construction  operations  were  monitore^d  by  a ae- 
^ tialed  instrumentation  program  including  p lexoi.ietei  s , 

I slope  inclinometers  and  settlement  plates  and  plug  ..  I i' 

iJ  zomotric  measurements  taken  during  construction  ui'i.  -i’. 

' that  little  pore  pressure  dissijjation  occurred.  Ff  r ‘ 

t! 

I case,  the  undrairu'd  slope  stability  analysis  is  u : , 

Prior  to  construct  ion  of  Test  Section  II 

I.; 

f*'  turbed  soil  borings  were  madt-  to  FI.  -120m-  ». 

1 I 

' were  located  at  the  levee  centerline  .ml  t* 

I of  f set s to  both  sides  of  the  cep. t i ; 1 . • . ,i  • 
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stations.  Continuous  samples  wore  taken  with  a 5-inch 
diameter  steel  tube  piston  type  sampler.  The  laboratory 
testing  program  included  visual  classification,  Attorborg 
limits;  unconfined  compression  tests  (UC) , unconsolidated 
undrained  (UUC)  and  isotropically  consolidated  undrained 
(CIU)  triaxial  compression  tests;  consolidated  drained  (CD) 
direct  shear  tests;  lab  vane  tests;  undrained  creep  tests; 
and  consolidation  tests.  The  undrained  shear  strength  was 
also  measured  in  the  field  with  a 2-inch  diameter  hand 
cranked  vane  test.  Test  holes  were  drilled  adjacent  to 
each  undisturbed  boring.  The  vane  strength  tests  were  per- 
formed at  5-foot  increments  of  depth.  Boring  locations  are 
shown  in  Figure  3-10. 

Due  to  the  proximity  of  the  adjacent  borrow  pit, 
floodwayside  stability  is  more  critical  than  landside  sta- 
bility. Therefore,  this  case  study  will  focus  on  the  flood- 
wayside  stability  of  the  main  embankment  at  TS  III.  A 
sharp  increase  in  strength  at  El.  -45msl  was  detected  vnth 
the  field  vane.  This  increase  was  found  to  be  the  result 
of  an  overconsolidated  layer  beginning  at  that  elevation 
(Foott  and  Ladd,  1973).  El.  -45msl  will  servo  as  a limiting 
depth  of  assumed  fciiluro  arcs.  Only  strength  measurements 
above  this  elevation  will  be  used. 

The  analyses  of  Test  .Section  III  is  subdivided  into 
two  i|rouf)S.  One  group  will  be  based  solely  on  field  vane 
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data,  the  other  on  only  unconfined  compression  data.  The 
reliability  of  the  main  embankment  of  Test  Section  III 
will  be  expressed  in  terms  of  safety  factors  and  failure 
probabilities. 

Figures  3-11  through  3-16  present  the  results  of 
strength  tests  performed  on  seimples  from  the  centerline, 
105ft  f loodwayside,  and  180ft  floodwayside  borings.  Fig- 
ures 3-11  through  3-13  show  the  average  results  of  field 
vane  testing;  Figures  3-14  through  3-16  show  the  average 
results  of  unconfined  compression  testing.  The  strength 
linos  presented  are  the  average  of  results  over  5ft  layers 
from  the  six  borings  at  each  respective  offset. 

Based  on  these  average  undrained  strengths,  the  sta- 
bility of  the  main  (floodwayside)  embankment  was  checked 
with  a circular  arc  method  of  analysis.  The  critical 
mean  safety  factors  for  field  vane  and  for  unconfined  com- 
pression strength  measurements  are  shown  in  Table  3-3. 

Also  shown  in  Table  3-3  are  the  critical  mean  safety  fac- 
tors corrected  for  bias  by  the  revised  Bjerrum  and  the  Wu 
correction  factors. 


CASE 

FS 

CF  (PI  = 70) 

FS  X CF 

FV 

2.08 

0. 69 

UC 

1.42 

1.10 

.STABILITY  OF  TS  III  MAIN  (FLOODWAYSIDE) 
EMBANKMENT  (Foott  and  Ladd,  1973) 

Table  3-3 
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From  the  results  of  conventional  stability  analyses, 
both  methods  of  measuring  strength  predict  satisfactory 
performance.  The  corrected  safety  factors  using  field  vane 
tests  and  unconfined  compression  tests  to  measure  strength 
were  within  fifteen  percent  of  each  other.  The  critical 
mean  safety  factor  based  on  the  unconfined  compression 
test,  when  corrected , was  higher  than  that  based  on  a field 
vane  tests. The  critical  failure  arcs  for  Test  Section  III 
are  shown  in  Figure  3-17  and  3-18. 

3.3.2  Probabalistic  Analysis 

The  corrected  safety  factors  presented  in  Table  3-3 
are  the  mean  safety  factors  of  the  main  (floodway side)  em- 
bankment using  two  different  methods  of  strength  deter- 
mination. The  failure  probabilities  of  this  embankment, 
calculated  with  the  different  strength  measurements,  are 
dependent  on  the  uncertainty  of  bias,  random  error  and  in- 
herent variability. 

BIAS  The  uncertainty  of  strength  bias  is  measured  by  the 
coefficient  of  variation  of  the  two  bias  correction  fac- 
tors. For  the  unconfined  compression  tests,  the  Wu  correc- 
tion factor  was  used.  This  correction  factor  has  a coeffi- 
cient of  variation  V , equal  to  0.30,  see  Figure  3-6. 

b 

The  clay  foundation  for  Test  Section  III  had  an  aver- 
age plasticity  index  equal  to  70  percent.  From  Figure  3-5, 
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the  revised  Bjerrum  correction  factor  for  this  typo  soil 
is  0.69.  The  coefficient  of  variation  due  to  this  bias  is 
equal  to  0.22  for  the  field  vane  test  results. 

RANDOM  TESTING  ERROR  The  spatial  coefficient  of  variation 
due  to  random  testing  error  is  calculated  by  equation 
(2.6-7b).  The  average  point  coefficient  of  variation  of 
strength,  Vs,  is  equal  to  0.32  and  0.28  as  measured  by  un- 
confined compression  tests  and  field  vane  tests,  respec- 
tively. The  mean  strength  from  unconfined  compression  tests 
reflect  the  average  of  239  tests;  the  mean  strength  from 
field  vane  tests  was  determined  from  178  tests. 

The  spatial  coefficients  for  variation,  Vg//n,  due 
to  random  testing  error  are  equal  to  0.02  for  both  methods 
of  testing.  This  source  of  uncertainty  is  small  relative 
to  the  uncertainty  of  the  bias  correction. 

INHERENT  SOIL  VARIABILITY  The  spatial  coefficient  of 
: variation  due  to  inherent  soil  variability  is  calculated 

by  equation  (2.6-16b).  The  critical  failure  arc  using 
' field  vane  strengths  was  approximated  with  three  perpondi- 

I cular  lines  as  shown  in  Figure  3-20.  The  lines  were  selec- 

I ted  to  preserve  both  the  depth  of  failure  and  the  length 

! of  the  failure  arc.  The  axial  length  of  the  critical  fail- 

< 

*•  ure  surface  was  set  equal  to  three  times  the  base  width, 

i 

Based  on  the  surface  area  of  the  simplified  failure 
plane  and  on  cor relat’ on  distances  selected  equal  to  5 
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feet  vertically  and  100  feet  horizontally,  the  number  of 
equivalent  soil  elements  is  determined  from  Figure  2-6. 

The  total  number  of  equivalent  elements  is  82.  For  this 
number  of  elements,  the  spatial  coefficient  of  variation, 
Vs/v'Ne,  due  to  soil  heterogeneity  is  equal  to  0.03. 

Through  similar  calculations,  the  spatial  coefficients  of 
variation  due  to  soil  heterogeneity  using  unconfined  com- 
pression strength  is  also  equal  to  0.03. 

A summary  of  the  data  for  Test  Section  III  is  shown 
in  Table  3-4.  Again,  the  uncertainty  of  the  bias  correc- 
tion is  the  main  source  of  uncertainty.  The  total  uncer- 
tainty in  the  predicted  average  shear  strength,  and  hence 
the  resisting  moment,  is  equal  to  0.22  and  0.30  for  field 
vane  testing  and  unconfined  compression  testing,  respec- 
tively. 

Based  on  the  above  uncertainties  and  the  safety  fac- 
tors shown  in  Table  3-3,  the  failure  probability  with  field 
vane  strength  measurements  is  calculated  as  9 percent 
(6  = 1.37)  using  equation  (2.4-1).  The  failure  probability 
with  unconfined  compression  tests  is  12  percent  (R  = 1.20). 
The  graphical  representation  of  those  failure  probabilities 
13  shown  in  Figure  3-21.  Also  shown  in  this  figure  are 
plots  of  failure  proliability  versus  safety  factor,  assum- 
ing that  the  uncertainty  in  the  resisting  moment  is  con- 
stant. 
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It  is  interestincj  to  note  that  the  safety  factor 
calculated  from  field  vane  data  is  lower  than  that  from 
unconfined  compression  data;  however,  the  probability  of 
failure  calculated  from  field  vane  data  is  lower  than  that 
from  unconfined  compression  data.  This  is  due  to  the  in- 
creased uncertainty  in  applying  the  Viu  bias  correction 
factor . 

To  investigate  the  importance  of  failure  arc  length 

on  the  calculated  failure  probability,  a smaller  failure 

arc  having  a higher  safety  factor  based  on  field  vane 

strengths  was  checked.  The  comparative  failure  arcs  are 

shown  in  Figure  3-18.  The  uncertainties  due  to  bias  and 

to  random  testing  error  are  the  same  for  both  failure  arcs. 

However,  the  spatial  uncertainty  due  to  soil  variability 

increases  from  0.03  to  0.04  { N decreases  from  82  to  54 

e 

due  to  the  decrease  in  the  failure  surface  size) . This 
change  in  uncertainty  has  a negligible  effect  on  the  total 
uncerta inty . 

The  failure  probability  for  the  smaller  arc  is  7 
percent  ( (<  = 1.51).  Although  the  critical  failure  arc 
bas<'d  on  minimum  safety  factor  does  not  necessarily  define 
the  arc  with  the  maximum  failure  probability,  it  generally 
does  for  the  above  examples  due  to  the  significance  of 
the  uncertainty  in  bias.  The  data  for  the  smaller  failure 
arc  are  shown  in  Table  3-4. 
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Test  Section  III  did  not  f^il.  It  did  undergo  large 
deformations  and  cracking  along  the  embankment  crown. 

These  deformations  indicate  that  the  embankment  was  on  the 
verge  of  failure.  This  nearness  to  failure  is  not  surpris- 
ing in  view  of  the  high  failure  probabilities  based  on 
both  methods  of  strength  measurements. 
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Chapter  4 


The  Inf luence  of  Embankment  Length 


on  Failure  Probabilit'i 


4.1  General . In  Chapter  2,  a basic  probability  model  for 
analyzing  the  reliability  of  an  embankment  constructed  on 
soft  saturated  clay  was  presented.  This  chapter  attempts 
to  assess  the  effects  of  embankm-^nt  length  on  the  computed 
failure  probability.  Before  discussing  the  effects  of 
length,  one  must  reflect  on  what  the  failure  probability 
computed  from  the  basic  model  represents. 

The  basic  model  computes  the  failure  probability,  pf, 
at  one  "point",  which  may  be  thought  of  as  the  center  of  a 
segment  (along  the  embankment  centerline) , whose  minimum 
length  is  selected  to  satisfy  the  assumption  of  "plane 
strain".  The  model  also  assumes  that  all  strength  measure- 
ments are  made  on  the  critical  failure  surface  for  the 

length,  L.  For  the  FRT  test  section,  this  minimum  length 
( * ) 

is  240  feet  , whereas  for  TS  III,  the  minimum  length  of 

( * ) 

embankment  is  810  feet 


If  the  actual  length  of  embankment  constructed.  La, 
were  equal  to  these  selected  lengths,  then  the  computed 
failure  probability  would  correctly  estimate  the  embankment 


These  minimum  lengths  wore  computed  in  Chapter  3.  By  mocJ(?l- 
inq  tht*  critical  failure  arc  with  a rectangular  surface,  tJie 
minimum  length  was  assumed  to  be  three  times  the  base  width 
of  the  simplified  failure  surface.  Those  minimum  lengths  are 
presented  in  Figures  3-7  and  J-i9. 
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failure  probability.  However,  if  the  actual  embankment  wore 
longer,  the  computed  failure  probability  would  tend  to  be 
unconservative  on  account  of  the  series  system  or  "weakest 
link"  concept  of  reliability. 

Two  approaches  to  handling  length  effects  are  pre- 
sented in  the  following  two  sections.  The  first  is  an  exten- 
sion of  the  basic  model  based  on  system  reliability  theory. 
The  second  is  a three-dimensional  model  developed  from  a 
"first  passage"  failure  criterion. 


4.2  Series  System  Failure  Criterion.  As  pointed  out  above, 
if  an  embankment  is  constructed  just  long  enough  for  one 
failure  to  occur,  then  the  basic  model  gives  a correct  esti- 
mate of  failure  probability. 

As  the  length  of  the  actual  embankment.  La,  exceeds 
the  minimum  length,  L,  the  estimated  failure  probability 
will  increase.  Considering  the  actual  embankment  length  as 
a multiple  of  the  minimum  length  (n  = La/L)  and  by  consider- 
ing the  adjacent  lengths  as  independent  events  having  the 
same  computed  failure  probability,  the  probability  of  a 
failure  occurring  somewhere  along  the  entire  length  can  be 
expressed  approximately  as  (Benjamin  and  Cornell,  1970).' 


or 


= 1 - (1  - pf) 


(if  n , •'<  1) 


(4.2-la) 
(4 . 2-lb) 
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where 


= the  probability  of  a failure  occurring 
at  some  location  along  the  entire  em- 
bankment length,  La. 

pf  = the  probability  of  a failure  occurring 
at  one  "point" , centered  at  the  minimum 
embankment  length,  L. 

n = the  number  of  minimum  lengths  contained 

( * ) 

within  the  actual  embankment  length 
The  implication  of  equation  (4.2-1)  can  be  shown  by 
reviewing  the  failure  probabilities  of  the  two  case  studies 
presented  in  Chapter  3.  Considering  the  calculations  based 
on  field  vane  testing,  the  Fore  River  Test  Section  had  a 
computed  failure  probability  of  0.4  percent  for  a minimum 
length  equal  to  240  feet.  Since  the  FRT  test  section  was 
itself  approximately  240  feet  square,  n equals  4 for  the 
entire  length  susceptible  to  failure.  The  probability  of  a 

failure  occurring  on  one  of  the  four  sides  of  the  test  em- 

(** ) 

bankment  is  calculated  as  1 . 6 percent  from  equation  (4.2-1). 


* 

This  model  is  based  on  the  assumption  that  the  average 
strengths  in  adjacent  segments  are  statistically  independent. 
Since  horizontal  correlation  distances  (assumed  equal  to  100 
feet  for  case  studies  in  Chapter  3)  are  usually  much  smaller 

than  the  minimum  lengths,  this  assumption  appears  justified. 

★ * 

The  clay  layer  under  two  of  the  sides  (South&West)  of  the 
test  section  is  not  as  deep  and  consequently  not  as  likely  to 
fail  as  the  North  and  East  slope.  This  effect  of  geometry  is 
not  considered  in  this  analysis.  It  is  assumed  that  the 
layer  thickness  is  constant  and  equal  to  30  feet. 
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For  the  Atchafalaya  Basin  Test  Section  III,  the  con- 
structed length  was  1500  feet,  while  the  estimated  minimum 
length  for  a "plane  strain"  failure  was  810  feet.  The  pro- 
bability of  a failure  occurring  along  the  minimum  length, 
based  on  field  vane  tests,  was  computed  as  9 percent.  This 
is  compared  to  a failure  probability  for  the  entire  length 


equa  1 

to  17  percent. 

The  results 

of  these  two 

case  studies 

are  shown  on  Table  4 

-1. 

Case 

Study 

Minimum 
Length, L 
(ft) 

Actual 
Length , La 
(ft) 

★ 

pf  for  L 
(%) 

★ 

P,for  La 
(%) 

1 

__  _ . _ . 

FRT 

240 

960 

0.4 

1.6 

TS  III 

810 

1500 

9 

i 

*Based  on  Field  Vane  Tests 

Failure  Probabilities  Including  Length 

Effect 

TABLE  4-1 


4.3  "First  Passage"  Failure  Criterion.  An  alternative  to 
the  above  method  of  introducing  the  parameter,  embankment 
length,  into  a probabilistic  analysis  of  stability  is  devel- 
oped from  a "first  passage"  failure  criterion. 

For  this  analysis,  the  standardized  safety  margin, 
is  now  thought  of  as  a stationary  Gaussian  random 
process : 
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^x)  = 


a 


M 


whore 


M,  = M , , - M , , 
(x)  r(x)  o(x) 


(4.3-1) 


and 

M . = a "moving  average"  of  strength  as  a 

r ^ X j 

function  of  position  along  the  embar.V.- 
ment  length. 

= a "moving  average"  of  the  load, 
assumed  constant. 


The  definition  of  U,  . is  illustrated  in  Figure  4-1.  The 

IX/ 

process  begins  and  ends  at  a distance,  L/2,  from  both 

ends  of  the  emban)cment.  Therefore,  at  x = 0,  the  actual 
emban)cment  length.  La,  is  equal  to  the  minimum  length  re- 
quired for  a "plane  strain"  failure  to  occur. 

The  probability  of  a failure  occurring  at  some  loca- 
tion along  a constructed  embankment  can  be  defined  as  the 
probability  that  the  standardized  safety  margin  lies  below 
the  reliability  index  at  any  point,  or  that  the  length  to 
first  passage  of  the  reliability  index  is  less  than  the 

length  of  the  process,  U,  . expressed  as: 

( X ; ^ 

^f  " ^ ^ ^(x)  f "PI  = P XI  (4.3-2) 

where 

Xp  = the  length  to  first  passage  of  the' 
barrier,  -B  (failure). 
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A widely  used  approximation  (Vanmarcke,  1972),  based 
on  the  assumption  that  crossings  of  a sufficiently  high 
barrier,  -B,  occur  independently  according  to  a Poisson 
process  with  rate,  v„  is; 

P 9 

P = (1  - pf)  exp  {-v_X}  (4.3-3a) 

s tJ 

Pg  = (1  - pf)  exp  exp(-6V2)}  (4.3-3b) 

where 

Pg  = the  probability  of  complete  survival  of  the 
embankment . 

pf  = the  probability  of  a failure  centered  at  a 
randomly-chosen  point. 

Vg  = the  mean  rate  of  uperossings  (i.e.,  crossings 

with  a positive  slope)  of  a level,  -6,  by  a 

stationary  Gaussian  process,  U,  . 

( X j • 

= the  moan  rate  of  uperossings  of  the  level 
" zero" . 

X = the  length  or  duration  of  the  process,  U.  . 

V ^ ^ » 

equal  to  La  - L. 

These  uperossing  rates  and  the  length,  X,  are  shown  in 
Figure  4-1. 

The  probability  of  a failure  occurring  somewhere 
along  the  embankment  length  is  equal  to: 

Pf  = 1 - Pg  (4.3-4a) 
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Pj  = 1 - {1  - pf)  exp  exp(-6^2)} 


(4. 3-4b) 


or  ~ pf  + v^X  exp(-6^/2),  if  sum  <<  1 (4.3-4c) 


By  inspecting  equation  (4.3-4b),  it  can  be  seen  that  as  X 
equals  zero  (the  actual  embankment  length.  La,  equals  the 
minimum  length,  L) , then 

Pf  = pf  (4.3-5) 

This  result  from  the  three-dimensional  model  is  consistent 
with  the  interpretation  of  the  result  of  the  basic  proba- 
bility model. 

The  major  difficulty  with  this  approach  to  length 

effects  lies  in  the  evaluation  of  v By  definition,  v is 

0 . 0 

a measure  of  how  rapidly  the  standardized  safety  margin 
fluctuates  with  position  along  the  embankment  length.  Con- 
sidering the  embankment  load  as  deterministic  and  constant, 
soil  strength  and  foundation  geometry  are  the  major  sources 
that  cause  the  standardized  safety  margin  to  fluctuate. 

Considering  only  the  soil  strength  as  variable,  what 
is  important  is  the  spatial  average  shear  strength  over  the 
critical  failure  surface  (three-dimensional).  This  spatial 
average  shear  strength  should  be  highly  correlated. 

Referring  back  to  Figure  4-1,  the  first  possible 
failur<-  surface  with  length,  L,  has  some  spatial  average 
•shf'.ir  .strength  acting  on  it.  This  average  shear  strength 
can  be  thought  of  as  centered  at  point  c^  Now  consider 
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another  possible  failure  surface  centered  at  some  small 


distance  from  c^  and  having  the  same  minimum  length,  L. 

The  latter  failure  surface  has  nearly  the  same 
spatial  average  shear  strength  as  the  former.  This  is  true 
since  a large  portion  of  the  two  failure  surfaces  are  common 
to  each.  It  is  this  overlapping  of  individual  failure  sur- 
faces when  computing  a "moving  average",  together  with  the 
large  horizontal  correlation  distances  estimated  for 
strength,  that  cause  the  spatial  average  shear  strength  to 
be  highly  correlated. 

The  mean  rate  of  zero  upcrossings  of  the  average 
process,  U,  . is  linked  to  the  correlation  distance  of  the 
spatial  average  shear  strength  (6<s>)  as  averaged  over 
three-dimensional  failure  surfaces.  If  the  coefficient  of 
correlation  of  this  spatial  average  shear  strength  can  be 
considered  to  have  the  form: 

p = g-x^/6s^  (4.3-6) 

then  the  rate  of  zero  upcrossings  becomes: 


Evaluating 
average  strength 
in  this  thesis, 
equation  (4.3-4b) 
Atchafalaya  Basin 


V = 


2/2 


(4.3-7) 


o 6<s> 

the  correlation  distance  of  the  spatial 
is  a difficult  task  and  is  not  undertaken 
However,  to  illustrate  the  application 

, different  values  of  v are  tried.  The 

o 

Test  Section  III  analysis  is  used  as  an 
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example.  The  probability  of  failure  of  a minimum  length  of 
this  embankment  (810  ft)  had  a computed  failure  probability 
of  failure  equal  to  9 percent  based  on  field  vane  testing. 

Figure  4-2  presents  plots  of  failure  probability 
versus  embankment  length.  Curve  1 is  based  on  the  three- 
dimensional  extension  presented  in  Section  4.2  (equation 
4.2-1).  Curve  2 is  computed  from  equation  (4.3-4b)  with 
set  to  yield  the  same  results  as  curve  1.  Curve  3 is 
also  computed  from  equation  (4.3-4b)  with  equal  to  1/L 
for  this  case. 

The  results  indicate  that  both  methods  of  predicting 
the  effect  of  embankment  length  on  failure  probability  can 
be  adjusted  to  yield  the  same  results.  From  curves  1 and  2 
on  Figure  4-2,  doubling  the  embankment  length  roughly 
doubles  the  failure  probability. 

Curve  3 is  based  on  a mean  rate  of  upcrossings  four 
times  greater  than  that  for  curve  2.  Using  this  greater 
value  of  upcrossinc  rate  generally  increases  the  failure 
probability  by  a factor  of  two  over  those  from  curves  2 and 
3. 

At  present,  it  cannot  b<,-  said  as  to  which  approach 
to  length  effect  is  better.  More  research  is  needed, 
especially  in  the  area  of  spatial  strength  variability.  Also, 
the  effect  of  foundation  geometry  (e.g.. 


natural  variability 


ot  depth  of  clay),  neqlcctcd  in  this  chapter,  may  be  a 
major  factor  influencing  failure  probability  versus 
length. 
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Chapter  5 Conclusions  and  Recommendations 


5 . 1 Genera  I 

A probabilistic  approach  to  analyzing  end-of-con- 
struction  embankment  stability  on  soft  saturated  clay  was 
presented  in  this  thesis.  The  model  is  proposed  as  a sup- 
plement to  the  current  design  approach  based  on  (J>  = 0 ana- 
lyses using  the  undrained  strength  of  the  foundation  clay. 
This  chapter  summarizes  the  conclusions  drawn  from  this 
study  and  also  makes  some  recommendations  for  future  re- 
search . 

5 . 2 Conclusions 

5.2.1  Need  to  Recognize  Uncertainty 

Chapter  one  presents  the  limitations  of  expressing 
the  reliability  of  an  earth  embankment  in  terms  of  a con- 
ventional safety  factor.  The  [principal  limitations  lie  i;i 
the  inability  of  the  current  design  approach  to  account  for 
the  uncertainty  inherent  in  the  approach  and  to  express 
the  risk  implied  in  design  in  quantitative  terms.  The  mag- 
nitude of  the  uncertainty  in  the  prediction  of  stability 
is  a critical  issue  in  geotechnical  engineering. 

It  is  the  role  of  the  geotechnical  engineer  to  reduce 
unk uowns  and  to  evaluate  and  assess  the  degree  of  uncertain- 
ty. The  probability  model  presents  a systematic  method  of 

qua  1 i ta t i vely  evaluating  uncertainty.  It  provides  the  en- 

-132- 


gineer  with  a mechanism  that  can  indicate  where  additional 
investigation  is  needed  to  min  imi  7,0  risk. 

The  general  probability  model  developed  in  this  the- 
sis recognizes  three  independent  sources  of  uncertainty  in 
undrained  strength  --  bias,  random  testing  error,  and  in- 
herent soil  variability.  By  evaluating  each  source  or  un- 
certainty independently,  the  engineer  is  able  to  see  where 
the  major  source  of  uncertainty  lies  and  concentr.ite  his 
efforts  to  minimze  it. 

Without  recognizing  and  treating  uncertainty,  the 
true  implication  of  safety  factors  are  unknown.  At  best,  the 
safety  factor  can  be  a tool  for  ranking  alternative  embank- 
ment designs  according  to  relative  reliability.  The  true 
meaning  of  the  safety  factor  varies  depending  on  many  fac- 
tors, such  as  how  s^  is  measured,  number  of  determinations, 
how  design  s is  selected  relative  to  average  s , and  me- 
thod  of  computing  factor  of  safety,  etc. 

Attempts  at  quantifying  the  uncertainty  involved  in 
analyzing  embankments  have  only  recently  been  made  (Barba- 
teu,  1972;  Yucemen  et.  al.,  1973;  Wu , 1974).  Correction 
factors  have  been  published  for  both  field  vane  tests 
(Bjcrrum,  1972)  and  unconfincd  compression  tests  (Wu,  1974) 
Those  mean  correction  factors  and  the  data  from  which  they 
wore  developed  allow  one  to  evaluate  the  uncertainty  asso- 
ciated with  applying  them. 
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A Study  is  currently  underway  at  M.I.T.  to  evaluate 
the  spatial  correlation  properties  of  undrained  shear 
strength.  Through  such  a study,  a better  understanding  of 
the  uncertainty  due  to  inherent  soil  variability  will  be 
available. 

Although  more  refinement  is,  indeed,  necessary,  the 
probability  model  is  a step  in  the  right  direction.  The 
probability  model  does  not  eliminate  any  of  the  steps  pre- 
sent in  the  current  design  approach.  But  is  does  attempt 
to  standardize  the  approach  and  it  introduces  another  set 
of  parameters  that  truly  enable  assessments  of  reliability 
of  analyses  --  those  that  measure  uncertainty. 

5.2.2  Safety  Factors  as  a Measure  of  Reliability 

Historically,  safety  factors  have  been  used  as  mea- 
sures of  reliability.  In  order  to  decide  how  good  safety 
factors  are  at  measuring  reliability,  it  is  necessary  to 
I reflect  on  the  decisions  that  influence  the  determination 

I of  a safety  factor  of  an  embankment  constructed  on  soft 

saturated  clay. 

A design  engineer  makes  decisions  that  directly  in- 
■ fluonce  the  calculated  safety  factor  of  an  embankment. 

• These  decisions  include  selecting  an  appropriate  method 

of  analysis,  planning  the  site  exploration,  laboratory  and 
field  testing  programs,  and  selectinq  the  "design"  strength 
for  the  analysis. 
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These  decisions  generally  do  not  follow  a standar- 


dized pattern.  For  a given  embankment  of  soft  clay,  an 
engineer  may  use  a wedge  or  a circular  arc  analysis;  2'-s- 
inch  or  5-inch  diameter  undisturbed  samples;  field  vane  or 
unconfined  compression  testing;  and  the  use  of  mean 
strengths  or  lower  than  mean  strengths  for  design.  Differ- 
ent combinations  of  the  above  selections  can  result  in  a 
very  wide  range  of  safety  factors  for  the  same  embankment. 
Consequently,  the  fact  that  one  embankment  has  a computed 
safety  factor  equal  to  1.5,  while  that  of  another  equals 
1.2,  may  have  little  real  significance  regarding  the  true 
"safety"  or  reliability. 

However,  safety  factors  can  be  used  to  reflect  rela- 
tive reliability,  at  least  for  the  same  soil  deposit,  if 
engineers  or  government  organizations  develop  standar- 
dized procedures  for  investigating  the  stability  of  embank- 
ments. By  building  upon  the  experience  of  the  engineer  and 
by  adjusting  the  design  procedure,  an  engineer  can  develop 
a design  procedure  that  can  limit  the  number  of  actual 
failures  to  an  acceptable  percentage.  This  usually  entails 
a trial  and  error  approach.  As  one  procedure  leads  to  too 
many  failures  or  too  conservative  of  a design,  one  or  a 
number  of  the  stops  are  changed.  This  can  ultimately  lead 
to  a combination  of  one  method  of  analysis,  a standardized 
field  investigation  program,  one  type  of  test  for  strength 
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determination,  and  one  means  of  selecting  a design  value 
of  strength.  Such  a combination  plus  a design  safety  fac- 
tor can  be  successfully  used  as  a semi-empirical  procedure 
for  investigating  stability  problems. 

By  standardizing  the  design  procedure  and  by  applying 
it  to  a particular  geographic  area,  the  uncertainty  in  the 
procedure  is  approximately  constant.  If  the  uncertainty 
(measured  in  terms  of  coefficients  of  variation  and  safety 
margin)  is  constant,  the  reliability  index,  and  hence  the 
failure  probability,  is  proportional  to  the  safety  factor. 
For  this  case,  the  safety  factor  can  be  used  as  an  indica- 
tor of  the  reliability  of  an  embankment. 

However,  there  are  limitations  to  such  semi-empir- 
ical  procedures.  Although  a design  safety  factor  is  used, 
the  true  safety  factor  of  an  embankment  remains  unknown. 
Biases  in  strength  measurements  and  the  method  of  analysis 
are  hidden  in  the  steps  of  the  procedure.  The  working  pro- 
cedure often  involves  an  cancellation  of  errors. 

If  another  method  of  determining  strength  were  used 
to  better  define  the  true  in  situ  strength,  the  balance  of 
the  procedure  can  be  destroyed  and  the  probability  of  fail- 
ure will  change  - even  though  the  design  safety  factor  is 
constant.  Also,  a work  ing  procedure  developed  for  a parti- 
cular geologic  deposit,  if  extrapolated  to  different  geo- 
logic areas,  may  prove  to  bo  ineffective. 
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Although  a working  factor  of  safety  procedure  does 
allow  a ranking  according  to  relative  reliability,  safety 
factors  cannot  be  combined  with  the  cost  associated  with 
an  embankment  to  optimize  design.  The  effect  on  the  fail- 
ure probability  by  raising  or  lowering  the  safety  factor 
cannot  be  evaluated  without  explicitly  accounting  for  the 
uncertainty  entering  into  the  analysis.  Although  increasing 
the  safety  factor  decreases  the  probability  of  failure, 
one  is  unable  to  answer  the  question  "By  how  much?" 

While  the  current  design  approach  has  often  proved 
to  be  an  effective  tool,  a systematic  consideration  of  un- 
certainty can  give  better  measure  of  true  risk.  With  such 
an  assessment  of  true  risk,  the  engineer  is  better  equipped 
to  optimize  the  design  and  judge  its  acceptability. 

5.2.3  The  Uncertainty  of  Bias 

The  two  case  studies  presented  in  this  thesis  indi- 
cate that  the  bias  inherent  in  field  vane  and  unconfined 
compression  testing  can  be  the  major  source  of  uncertainty 
in  design.  Depending  on  soil  type  and  depth  and  type  of 
tost,  meaasured  strength  can  be  equal  to,  less  than,  or 
greater  than  the  true  in  situ  undrained  strength.  Bias  is 
a [)orsi!3tant  source  of  uncertainty  regardless  of  the  number 
of  1 nd('p«'ndont  strength  tf'St  s performed. 

Comparing  the  uncertainty  of  the  two  correction 
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factors  used  for  the  examples  presented  in  this  thesis, 
the  correction  factor  for  unconfined  compression  tests 
showed  the  larger  scatter  about  its  mean.  This  would  tend 
to  indicate  that  the  field  vane  test  is  a better  measure 
of  strength  trends  in  space  than  the  unconfined  compression 
test . 

The  reason  for  field  vanes  being  better  measures  of 
strength  is  that  the  field  vane  test  introduces  less  (and 
generally  a constant  amount  of)  sample  disturbance.  Vary- 
ing degrees  of  sample  disturbance  are  exhibited  by  uncon- 
fined compression  tests  by  different  size  and  type  of  sam- 
pling tubes,  different  field  and  lab  handling  techniques, 
different  type  soils  tested.  It  also  tends  to  vary  with 
depth.  The  field  vane  tests  minimizes  sample  disturbance 
by  performing  the  test  in  situ.  It  also  generally  involves 
a more  standardized  test  procedure. 

There  appears  to  be  two  different  ways  to  reduce  the 
uncertainty  due  to  test  bias.  The  first  involves  restric- 
ting the  source  of  data  for  the  correction  factor.  The 
second  way  is  to  use  a more  rational  method  of  measuring 
strength . 

By  limiting  the  source  of  data  to  a particular  geo- 
logic deposit,  the  effect  of  soil  type  on  the  correction 
factor  Will  be  approximately  constant.  By  standardizing 
both  the  equipment  and  the  method  of  testing  (i.e.,  strain 
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rate,  depth  of  penetration  of  the  vane,  size  of  the  vane) , 
the  bias  due  to  these  sources  of  error  will  also  be  ap- 
proximately constant. 

The  effect  of  restricting  the  method  of  testing  and 
source  of  data  for  the  correction  factor  was  illustrated 
in  Chapter  1.  By  limiting  the  source  of  case  studies  to 
embankment  failures  in  Boston  Blue  Clay  and  by  performing 
the  field  vane  test  in  a standardized  manner,  scatter  in 
a corrrection  factor  developed  from  this  information 
should  be  smaller  than  that  for  the  revised  Bjerrum  fac- 
tor (developed  from  a wide  variety  of  case  studies) . 

Of  course,  one  problem  with  developing  local  correc- 
tion factors  is  that  it  takes  time  to  accumulate  a suffi- 
cient amount  of  data.  Performing  strength  tests  at  sites 
of  embankment  failures  after  the  failures  have  occurred 
would  hasten  the  acquisition  of  local  data. 

The  second  way  to  reduce  the  uncertainty  due  to  bias 
is  to  improve  on  the  conventional  method  of  measuring  soil 
strength.  Such  a method  has  been  proposed  by  Ladd  and  Foott 
(1973)  and  has  shown  to  be  more  effective  in  predicting  un- 
biased strength  for  uniform  soil  deposits.  The  Stress  His- 
tory and  Normalized  Soil  Engineering  Properties  (SHANSEP) 
method  is  based  on  the  principle  that  stress-stra  in- 
strength  properties  of  clay  ar(>  unicjuely  related  to  the 
overconsolidation  ratio  (OCR).  By  reconsolidating  a sample 


in  the  laboratory  to  the  same  stress  ratio  as  in  the  field, 
one  can  obtain  a sample  having  normalized  properties  simi- 
lar to  those  in  situ.  By  shearing  the  reconsolidated  sam- 
ple under  the  same  stress  conditions  expected  in  the  field 
(plane  strain  active,  direct  shear,  plane  strain  passive), 
normalized  strength  parameters  (undrained  strength  divided 
by  vertical  effective  stress)  can  be  determined  for  a soil 
deposit . 

The  measurement  of  strength  by  the  SHANSEP  method 
does  not  involve  a direct  measurement  of  in  situ  strength 
at  a point.  Rather,  normalized  soil  properties  are  ob- 
tained for  an  entire  soil  deposit.  The  evaluation  of  the 
uncertainty  of  strength  determined  by  SHANSEP  requires  an 
evaluation  of  the  uncertainty  in  the  steps  followed  to 
measure  normalized  soil  properties  (total  weight,  pore 
pressures,  coefficient  of  lateral  earth  pressure,  maximum 
past  pressure,  laboratory  strength  testing) . Preliminary 
studies  in  this  area  indicate  that  uncertainty  in  stress 
history  (especially  maximum  past  pressure)  is  the  dominant 
source  of  uncertainty  for  strength  measured  by  SHANSEP. 

S . 3 Recommendations 

S . 3 . 1 Othe r Failure  Modes 

The  prol>ability  model  developed  in  this  thesis  as- 
sumes only  one  mode  of  embankment  failure,  an  undrained 


-140- 


shear  failure.  It  neglects  other  possible  inodes  of  failure, 
such  as  excessive  consolidation  deformations,  deformations 
due  to  undrained  creep,  piping,  and  overtopping  or  erosion. 

If  other  modes  of  failure  are  to  be  considered,  an 
assessment  of  the  probability  of  failure  for  each  indivi- 
dual made  has  to  be  made.  Once  the  individual  failure  pro- 
babilities are  assessed,  they  can  be  combined  using  basic 
probability  to  give  a failure  probability  resulting  from 
the  system  of  failure  modes. 

More  study  on  a probabal istic  basis  is  needed  on 
other  modes  of  failure  of  earth  embankments.  Without  a full 
understanding  of  these  other  modes,  th'e  actual  failure 
probability  is  indeterminate.  This  area  of  study  would 
also  involve  introducing  another  variable,  time,  into  the 
problem. 

5.3.2  Spatial  Variability 

The  uncertainty  due  to  inherent  soil  variability  was 
shown  to  be  proportional  to  the  point  coefficient  of  varia- 
tion of  in  situ  strength  and  inversely  proportional  to  the 

(!qui valent  number  of  independent  soil  elements,  N . In  order 

e 

to  ('valuate  , the  correlation  of  in  situ  strength  has  to 
be  measured.  Determining  the  correlation  of  in  situ  strength 
would  generally  require  a close,  well  tested  grid  of  borings 

A study  is  currently  underway  at  M.I.T.  to  evaluate 
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strength  correlation  distances  for  various  soil  deposits. 
More  research  is  needed  in  this  area  of  uncertainty  in 
order  to  better  understand  the  true  importance  of  inherent 
soil  variability. 

5.3.3  Other  Measures  of  Strength 

As  previously  stated,  the  Stress  History  and  Nor- 
malized Soil  Engineering  Properties  approach  has  been  ef- 
fectively used  to  measure  undrained  strength.  A study  of 
the  sources  of  uncertainty  in  the  approach  and  a method 
for  measuring  the  uncertainty  should  be  made.  This  would 
allow  for  the  results  from  SHANSEP  to  be  applied  in  proba- 
balistic  studies  of  embankment  reliability.  Also,  more 
studies  of  other  embankment  failures,  as  well  as  of  satis- 
factorily performing  embankments,  should  be  made.  This  will 
provide  a better  feel  for  the  results  of  probabalistic 
analyses . 
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NOTATION 


! 


I 


! 


cc 

CF 

COV 

FS 

M 

M 


(x,y) 


M 


o (x) 


M 


R 

M..  / V 
R (x) 

Pf 


S 

S 


(i) 

(neas) 

(i) 

(true) 


total  costs  of  construction 

consequences  of  failure  expressed  as  costs 

covariance  of  x and  y 

safety  factor 

safety  margin 

overturning  moment 

a "moving  average"  of  the  overturning  moment 
resisting  moment 

a "moving  average"  of  the  resisting  moment 

probability  of  failure  at  one  "point" 

probability  of  a failure  occurring  at  some 

location  along  the  actual  embankment  length 

measured  strength  at  point  i 
true  strength  at  point  i 


S 

u 

TEC 


U 


undrained  shear  strength 
total  expected  costs 

standardized  safety  margin  at  a "point" 

standardized  safety  margin  as  a stationary 
Gaussian  random  process. 

Variance  of  x 

coefficient  of  variation  of  x 
length  or  duration  of  the  process, 
moan  of  x 
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Xg  length  to  first  passage  of  the  barrier,  -6 

6 reliability  index 

6s  correlation  distance  of  strength 


random  testing  error  at  point  i 
bias  or  systematic  error  at  point  i 
Bjerrum  correction  factor 
variance  of  x 
standard  deviation  of  x 
Wu  correction  factor 

mean  rate  of  upcrossings  of  a level, 

mean  rate  of  upcrossings  of  the  level,  zero 
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